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Entanglement swapping is a process by which two initially independent quantum systems can become entan- 
gled and generate nonlocal correlations. To characterize such correlations, we compare them to those predicted 
by bilocal models, where systems that are initially independent are described by uncorrelated states. We ex- 
tend in this paper the analysis of bilocal correlations initiated in [Phys. Rev. Lett. 104, 170401 (2010)]. In 
particular, we derive new Bell-type inequalities based on the bilocality assumption in different scenarios, we 
study their possible quantum violations, and analyze their resistance to experimental imperfections. The bilo- 
cality assumption, being stronger than Bell's standard local causality assumption, lowers the requirements for 
the demonstration of quantumness in entanglement swapping experiments. 



I. INTRODUCTION 

The study of correlations between the outcomes of mea- 
surements performed on several quantum systems has led to 
remarkable progresses, both on fundamental aspects of quan- 
tum theory and on potential applications in quantum informa- 
tion technologies. Particularly intriguing are the cases where 
the various quantum systems are all at a distance from each 
other. From a fundamental point of view this situation led to 
the discovery of quantum nonlocality, that is, of the existence 
of correlations that cannot be described by a locally causal 
model (lfl. From an applied point of view, these studies led, 
quite recently, to the understanding of the power of nonlocal 
correlations for quantum information processing, in particu- 
lar for reducing communication complexity 12J], for Quantum 
Key Distribution (QKD) Q 0], private randomness genera- 
tion IHIH], or device independent entanglement witnesses 0]. 
Interestingly, in such examples nonlocal correlations can be 
exploited directly, in a device independent manner, indepen- 
dently of the Hilbert space machinery of the quantum theory. 
This applied side led, in turn, to a better understanding of 
some fundamental aspects of quantum theory and quantum 
information, such as for instance hidden assumptions in the 
abstract security analyses of QKD Hfl. 

The usual starting point in such works on nonlocality are 
limitations — such as, e.g., Bell inequalities |9[ — on the pos- 
sible correlations between the measurement results on distant 
systems, following from the principle of local causality. For- 
mally, the different systems measured in the experiment are 
considered to be all in an initial joint "hidden'fl state A, where 
A is arbitrary and could even describe the state of the entire 
universe prior to the measurement choices. The measurement 
outcome of any particular system can depend arbitrarily on the 
global state A and on the type of measurement performed on 
that system, but not on the measurements performed on dis- 
tant systems. This last condition is Bell's local causality as- 
sumption 0]] (or Bell's "locality assumption", simply), which 



We keep the terminology "hidden" to describe the states A for historical 
reasons QJ. Note that these states A need not actually be hidden, i.e., inac- 
cessible to the observer. 



implies, e.g. in the case of three parties, that the measurement 
outcome probabilities can be written as 

P(a,b,c\x,y,z) = JdXp(X)P{a\x,X)P{b\y,X)P{c\z,X), (1) 

where x, y, z denote the measurement settings ("inputs") cho- 
sen by the three parties, a, b, c denote the corresponding mea- 
surement outcomes ("outputs"), where p(X) is a probability 
distribution over the set of all possible joint hidden states A, 
and where it has implicitely been assumed that the measure- 
ment choices x, y, z are independent of A. 

Nowadays, fast progress towards advanced demonstrations 
of quantum communication networks, involvingquantum re- 
peaters n 1 Oh based on entanglement swappings 111 ill and quan- 
tum memories lfl2ll . are underway in many labs around the 
world. In these future quantum networks, several indepen- 
dent sources of entangled qubit pairs will distribute entangle- 
ment to partners who will then connect their neighbours by 
performing joint measurements on two (or more) qubits, each 
entangled with one neighbouring qubit, as illustrated for the 
simple case of three partners in Fig.Q] Such experiments have 
an interesting feature that has so far received little attention in 
previous works on nonlocality: the multipartite correlations 
between the measurement results at each site do not originate 
from a single multipartite entangled state, but from a series 
of bipartite entangled states that are initially independent and 
uncorrelated from each other; i.e., there is not a unique inital 
joint state (the analogue of A in a locally causal model) that is 
responsible for the observed correlations, but these are instead 
created from smaller systems through joint measurements. 

To understand and characterize the nonlocal properties ex- 
hibited in such experiments, it is natural to compare them to 
models where independent systems are characterized by dif- 
ferent, uncorrelated hidden states A. In the case, e.g., of the 
experiment of FigureQ] one would thus replace Bell's locality 
condition (Q~|i by 

P(a,b,c\x,y,z) = JJ dAidA 2 pi(Ai) p 2 (A 2 ) 

P{a\x : Xx)P{b\vMM)P{c\z,X % \ (2) 

where Ai characterizes the joint state of the systems produced 
by the source S\ and A 2 for the source 5* 2 ; see Figure [2] Of 
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FIG. 1. Typical entanglement swapping scenario where three parties, 
Alice, Bob and Charlie, share two sources Si and 52 that each emit 
independent pairs of particles in some quantum states gi and Q2- Bob 
performs a joint measurement y on the two particles he receives from 
each source and obtains an output b. Depending on Bob's outcome, 
Alice and Charlie's systems end up in one out of different possible 
entangled states. Alice and Charlie apply some measurements x and 
z on their particle and obtain outputs a and c. Such an experiment is 
characterized by a joint probability distribution P(a,b,c\x,y,z). 
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FIG. 2. The natural counterpart of the entanglement swapping sce- 
nario of Fig.Q]in terms of a locally causal model with two indepen- 
dent sources of hidden states: the systems produced by the source 5*i 
are characterized by hidden states Ai, while those from the source 5*2 
are characterized by hidden states A2. The two sources are assumed 
to be independent, hence the joint distribution p(\i, A2) of hidden 
states has the product form p(\i, A2) = p\ {\\)pi(\2), as in eq. l(2j. 



course, one can never exclude on pure logical grounds that 
systems that appear independent to us, such as pairs of parti- 
cles produced by different sources, are not in fact correlated 
in some hidden way. But, quoting Bell, "this way of arranging 
quantum mechanical correlations would be even more mind- 
boggling than one in which causal chains go faster than light. 
Apparently separate parts of the world would be deeply and 
conspiratorially entangled" [j]]. 

Motivated by the earlier works lfl3l [Till , the study of cor- 
relations between the results of measurements performed in 
quantum networks was initiated in a recent letter lfl5tl from 
the point of view just introduced. This leads to interesting new 
scenarios: because of the assumption that independent sources 
are characterized by different and independent A's, this low- 
ers the requirements on experiments for the demonstration of 
quantumness of a network and, to start with the simplest case, 
the demonstration of quantumness of an entanglement swap- 
ping process. Such studies may lead to new applications, in 
the spirit of device independent quantum information process- 
ing M- 

The general approach considered here should also con- 
tribute to the characterization of the nonlocal properties as- 
sociated to joint measurements, a question that has received 
little attention in traditional works on nonlocality (see how- 
ever Refs. lfl6l[l7ll ). Indeed, the joint measurements needed 
to connect neighbouring qubits in quantum networks are nec- 



essarily entangling, i.e., entanglement between remote parties 
appears through joint measurements and not from a joint state 
of distant systems — formally, the eigenvectors of the opera- 
tors that describe joint measurement are entangled. Recall 
that a joint measurement, the so-called Bell state analyzer, is 
also at the core of the celebrated quantum teleportation pro- 
tocol 0- Entanglement, the characteristic property of quan- 
tum mechanics in Schrodinger's words lfl9ll . thus plays a dual 
role, once allowing joint states of several systems and once 
allowing joint measurements. 

Before discussing in more details the results presented in 
this paper, let us finally stress how natural the assumption of 
independent A's is. Actually, an equivalent assumption is al- 
ready implicit in all standard tests of Bell inequalities. In- 
deed, in such tests one needs to assume that the measure- 
ment settings are random and independent of the entangle- 
ment source |H]; this is achieved, e.g., by having a local 
Quantum Random Number Generator (QRNG) determining 
the random settings [20]. But this makes sense only if one 
assumes that the sources in the QRNG are independent of the 
entanglement source and that they are characterized by inde- 
pendent and uncorrected A's. Consequently, our assumption 
is actually not new, but merely formalizes a usually tacit as- 
sumption and extends its scope to more advanced topologies 
of quantum networks. 



Structure of the paper 

In this article, we develop and formalize the approach in- 
troduced above for the simplest case: that of three partners on 
a line, with two independent sources of entangled qubits, as 
illustrated in Fig.Q] Following 11511 . we call bilocal the corre- 
lations that can be described as resulting from measurements 
of two independent hidden states Ai and A2, each produced by 
one of the two sources. Conversely, a correlation that cannot 
be described in such a way is called non-bilocal. 

We first introduce more formally the concept of bilocal- 
ity in Section [II] and show how the bilocality assumption can 
be tested. A first approach, developped in subsection III C| is 
to look for explicit bilocal decompositions; we introduce ef- 
ficient representations for (bi-)local models, which help the 
search for explicit decompositions, and where the bilocality 
assumption takes a very simple form. Another approach is to 
test Bell-like inequalities, which are satisfied by bilocal cor- 
relations but can be violated by non-bilocal correlations; in 
subsection III Dl we derive such (nonlinear) bilocal inequali- 
ties (eqs. d20l ), d23l ) and d27l)), which typically take the form 



\J\ < 1, 



(3) 



where / and J are linear combinations of probabilities 
P(a, b, c\x, y, z) (see eqs. (HMEJ, (ErES and ^3®). 

Section [III] is devoted to the study of how the correla- 
tions produced in quantum entanglement swapping experi- 
ments violate our bilocal inequalities. We analyse the situa- 
tions where Bob performs a complete Bell state measurement 
in subsection MI AI and where he performs partial Bell state 
measurements that allows him to distinguish different pairs 
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of Bell states (subsection HUB 11 1, or to distinguish two Bell 
states, while the other two give the same outcome (subsec- 
tion lHIB 2k In all cases we study the resistance to white noise 
of the quantum violations, and find that the required visibili- 
ties for demonstrating non-bilocality are significantly lower 
than for demonstrating Bell nonlocality. The resistance to de- 
tection inefficiencies is also analysed in some simple cases for 
the complete Bell state measurement in subsection lHICI 

In Section [IV] we address further issues on quantum non- 
bilocality. We study a trade-off between the resistance to 
noise of nonlocality and non-bilocality for quantum correla- 
tions (subsection II V At . and show that the two are not neces- 
sarily correlated but that the maximization of one is made at 
the expense of the other. We then investigate possible viola- 
tions of the bilocality assumption using non-maximally entan- 
gled states (subsection llVBt . We also address the question of 
classically simulating noisy entanglement swapping correla- 
tions, and introduce two protocols (with and without commu- 
nication) for that in subsection llV CI 

Finally, in Section [V] we come back to our justification for 
the assumption of independent sources, and to the idea that 
it is actually already implicitly used in standard Bell experi- 
ments. We illustrate this claim by showing that the assump- 
tion of local causality with independent sources in a bipartite 
Bell test is equivalent to an assumption of trilocality in a four- 
partite experiment. 



II. CHARACTERIZING BILOCAL CORRELATIONS 
A. The bilocality assumption 

We consider the scenario depicted in Figure [2] with three 
parties sharing two sources of independent hidden states; this 
is the simplest case where our assumption of independent 
sources of hidden states makes sense. 

In such a tripartite scenario, Bell's locality assumption |Q]] 
reads (as we recalled in eq. (Q3) 

P(a,b,c\x,y,z) = JdXp(X)P(a\x,X)P{b\y,X)P(c\z,X). (4) 

Here given the state A, the outputs a, b and c of the three par- 
ties, for inputs x, y and z, are determined respectively by the 
local distributions P(a\x, A), P(b\y, A) and P(c\z,X). The 
hidden states A follow the distribution p(X), normalized such 
that JdAp(A) = 1. 

If we now assume that the response of the three parties de- 
pend only on the states Ai or A2 characterizing the systems 
that they receive from the sources Si or S2, respectively, we 
write: 

P(a,b,c\x,y,z) = JJ dXidX 2 p{X 1 , A 2 ) 

P(a|x,Ai)P(6|y,A 1 ,A 2 )P(c|z,A 2 ), (5) 

Note that without making any further assumption, equation (O 
is equivalent to ©; indeed, p(X%, A2) could be different from 
zero only when Ai = A2 = A to recover ©. 



Now we introduce our crucial assumption: since the two 
quantum sources Si and S2 are supposed to be independent, 
we assume that this property carries over to the local model, 
and therefore the distribution of the hidden states Ai and A2 
should factorize, in the form 

p(Xi, A 2 ) = Pi(Ai)/9 2 (A 2 ) ■ (6) 

Together with (O, this defines our assumption of bilocality, as 
already expressed in equation (|2j. The hidden states Ai and 
A 2 now follow independent distributions pi(Ai) and p 2 (A 2 ), 
such that JdAi pi(Ai) = JdA 2 p 2 (A 2 ) = 1. 

Note that, as in the standard case of Bell locality, no re- 
strictions are made on the sets on which Ai, A 2 are distributed 
(apart from the fact that they must be measurable). Expanding 
the results of lfl5ll . we show in the subsections III CI and III Dl 
below that for finite numbers of inputs and outputs, eqs. © 
and (O lead nonetheless to implementable tests of bilocality, 
by looking for explicit bilocal decompositions or by testing 
(non-linear) Bell-like inequalities. Before that, let us briefly 
mention some general properties of the set of bilocal correla- 
tions. 



B. Topology of the bilocal set 

Any bilocal correlation is by construction also local; the set 
of bilocal correlations ("bilocal set", B) is therefore included 
in the set of local correlations ("local set", £): B C C. 

It is well known, and clear from the definition Q, that the 
local set is convex. On the other hand, because of the non- 
linear constraint (O, a mixture of bilocal correlations is not 
necessarily bilocal: the bilocal set is not convex. Determin- 
istic local correlations are bilocal; since they are the extremal 
points of the local set, C is actually the convex hull of B. 

One can further show (see Appendix[A} that the bilocal set 
is connected, and that its restriction to subspaces where the 
marginal probability distribution of Alice (or Charlie) is fixed 
is star-convex; star-convexity does not however hold for the 
whole bilocal set. 



C. Explicit bilocal decompositions 

From now on we consider scenarios with finite numbers of 
possible inputs and outputs. 

In that case, the local set C forms a convex poly tope l2lll . 
The description of the local polytope as the convex hull of 
a finite set of extremal points — corresponding to determinis- 
tic local correlations — allows one to use efficient numerical 
approaches based on linear programming to determine if a 
correlation P is local. Alternatively, the local polytope can 
be described in terms of its facets — corresponding to (possi- 
bly trivial) Bell inequalities — which can be enumerated algo- 
rithmically for a small enough number of inputs and outputs. 
Hence, in order to determine whether a correlation P is local 
or not, one can either solve a linear programming problem, or 
check that all Bell inequalities are satisfied. 
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Because the bilocal set is not convex, one cannot use stan- 
dard Bell inequalities to distinguish the bilocal and the non- 
bilocal correlations. The bilocal set is significantly more diffi- 
cult to characterize than the local set. Still, similar approaches 
can be used: one can describe the question whether a given 
correlation P is bilocal or not as a non-convex feasibility 
problem, or alternatively (as we will see in the next subsec- 
tion lll Dl i. one can derive non-linear inequalities that are satis- 
fied by any bilocal point. However, we do not have a system- 
atic practical approach to obtain such "bilocal inequalities"; 
and even in the simplest case we consider, we do not have a 
complete set of inequalities that would be sufficient to define 
the set B. 

We start here by giving alternative formulations for the bilo- 
cality assumption (00, which will prove more handy to use 
for practical purposes, in particular when looking for explicit 
bilocal decompositions. 



1. Bilocal decompositions onto deterministic correlations, 
with weights 

Consider a local correlation P, written in the form (0. 
It is well known l22ll that Alice's local response function 
P(a\x, Ai) can (without any loss of generality) be taken to be 
deterministic, i.e., such that it assigns a unique measurement 
output a to every input x: any randomness used locally by Al- 
ice can indeed always be thought of as being included in the 
shared random variable Ai. For a finite number of possible 
measurement inputs and outputs, there is a finite number of 
such deterministic strategies corresponding to an assignment 
of an output a x to each of Alice's N possible inputs x. We la- 
bel each of these strategies with the string a = a± . . . a/v and 
denote the corresponding response function P & (a\x) = S a , ax 
(with 6 m .n = 1 if m = n, 8 mt „ = otherwise). Similarly, the 
response functions P(b\y, Ai, A2) and P(c\z, A2) can also be 
taken deterministic; we label the associated strategies f3 and 7 
and the corresponding response functions are Ps(b\y) = 5b, p y 
and Pj(c\z) = <5 C|7a . 

Integrating over the set of all pairs (Ai, A2) that spec- 
ify the strategies a, f3, and 7 for Alice, Bob and Charlie resp., 
we can write (0 as 

P(a,b,c\x,y,z) = ^ q^ P a {a\x)Pp{b\y)P^{c\z) (7) 

with q & *^ = jJ Al2 dAidA 2 p(Ai, A 2 ) > and 

a/37 

£a/37 = 1- Eq. (0 corresponds to the well-known de- 
composition of local correlations as a convex sum of deter- 
ministic strategies, where the weights g a/ j 7 can be understood 
as the probabilities assigned by the source to the strategies a, 
ft and 7. 

Since a is specified here by Ai and 7 is specified by A2, 
then (with obvious notations) U^A^V = = A g x A?, 



and 

l&i = J2pQa^ = UaIxA? dAidA 2 p(Ai,A 2 ), (8) 

te = S 7 fc 7 = JJ \i &xA 2 dAidA 2 p(Ai, A 2 ), (9) 

<?7 = Eate 7 = UaixA? dAidA 2 p(Ai, A 2 ), (10) 

where A 1 = U fi Ag and A 2 ~ U 7 A? are the state spaces of 
the variables Ai and A 2 . 

Let us now assume that P is bilocal. One can see from ([8]- 
|ToT > that the independence condition (0 implies that 

for all a, 7, q & ^ = q s <? 7 ■ (11) 

The interpretation is clear: the strategies a and 7 being deter- 
mined by two independent sources, their probabilities should 
be independent. 

Conversely, any correlation P(a,b,c\x,y,z) satisfying (0 
and ( fTTT i can be written in the form (0. Indeed, since 
937 = q&q^, we can write q & ^ = q a q^q^. In- 
serting this expression in (0 and defining P a ,^{b\y) = 
T,p<lp\a7 P p( b \y)> we then find that P(a,b,c\x,y,z) = 
Ea^ P & (a\x)P & ^{b\y)P^{c\z), which is clearly of the 
form (0. We thus conclude that a tripartite correlation is bilo- 
cal if and only if it admits the decomposition (0 with the re- 
striction ( fTTT i. 

Such a description of bilocal correlations is easier to deal 
with than the defining assumptions (0 and (0, as it involves 
only a finite number of coefficients q s p^. One can thus now 
determine whether a correlation P is bilocal by searching such 
weights <7 g) g 7 > 0, with the linear constraints that they must 
reproduce the correlation P as in (0, and the quadratic con- 
straints of bilocality ( fTTT i. 

An even more compact representation of bilocal decompo- 
sitions can however be given in the following way. 

2. Decompositions in terms of "correlators" ejjj. 

For a given local (or bilocal) correlation P, the decomposi- 
tion (0 is in general not unique; Eq. only imposes a lim- 
ited number of constraints on the weights g a( g 7 , which are not 
enough to fix all of them. When dealing with a local decom- 
position, it is convenient to use a parametrization that clearly 
separates the parameters that are fixed, and those that are in- 
ternal degrees of freedom of the local model. A nice way to 
do it is to transform the weights q s ^ into "correlators" e^.. 

Since we will use this approach in the context of bilocal- 
ity, we present it here in the 3-partite case; note however that 
this representation can be useful in more general studies of 
(non-)locality — not only of bilocality — and it can easily be 
generalised to any A^-partite case. Also, for simplicity we 
consider here a scenario with binary inputs and outputs; the 
generalization to other scenarios can be cumbersome but is 
rather straightforward (see subsection HID 31 and Appendix [Bl 
for instance for the case where Bob has only 1 possible input, 
and 4 or 3 outputs). 

For binary inputs and outputs, Alice, Bob and Charlie's 
strategies a,j3, 7 simply contain two bits: a = aoa±, etc. 
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Let us then define, for a particular local decomposition of P 
in terms of weights q a p~, and for i = ioi\,j = j'oii an d 
k = k ki £ {00,01, 10, 11}, the coefficients 

where a ■ i = a i a + ct^ii, etc. The set of coefficients {e^} 
is equivalent to the set of weights {q & p^} (eq. ( fT2b is actu- 
ally a discrete Fourier transformation), and ( TT2T > can easily be 
inverted to obtain 

Thus, both representations (in terms of coefficients q a ^ or 
ejjj.) can be used to unambiguously define the local decom- 
position. 

We call the coefficients e^g "correlators", for the following 
reason: defining for a strategy a(3j, A x = (— I)"*, B y = 
(-l fy, and C z = (-1) 7 % we have 

m = ( Aq° A'j 1 B ° B^ 1 Cq° C^* 1 ) {g a(5 _ } ■ (14) 

where the average value is computed with the weights q s p^ 
(note that for a deterministic strategy, all values A x , B y and 
C z are defined simultaneously). 

Now, when i ^ 1, j ^ 1 and fc 7^ 1 (with the notation 
1 = 11), at most one term A x , B y and C z per party appears 
non-trivially in the average above. Hence, this average value 
can be obtained directly from the correlation P, and the cor- 
responding correlator^] ejjj, are therefore fixed by P: 

for i 7^ l,j 7^ 1, and k 7^ 1, 

eg S = {A«A*B*>B?C*>C*) P (15) 

with now A a = (— l) ax , where a x £ {0, 1} is Alice's output 
for the input x, and similarly for the other two parties. The av- 
erage is now computed from the correlation P. For instance, 
one gets ei ,io,io = (A BqCq) p , and e i,oo,oo = (Ai)p. 
Note that in particular, with = 00, one has egoo = 1 by 
normalization. 

On the other hand, when i = 1, j = lorfc = 1, the 
average value in ( TPfl ) cannot be obtained from P, since the 
measurement results A$ and A\, Bq and B\, Co and G\ are 
incompatible; the corresponding correlators are internal de- 
grees of freedom of the local model, only constrained by the 
nonnegativity of q & ^, i.e., 

for all a, p, 7, g (_i)*W+7-S ^ > . (16) 



2 When expressing correlations, it is often more convenient to consider ±1- 
valued outputs. Throughout the paper, we'll use lowercase variables for bit 
values 0/1, and uppercase variables for the corresponding bit values ±1, 
as for instance in = (— l) ax . 

3 More generally, for more than two inputs per party, the fixed correlators 
are those where the corresponding product averaged in H4\ only contains 
terms that refer to at most one input per party. 



Coming back to the bilocality constraint: in the correlators 
representation, one can easily check that the condition (fTTT i 
translates into the following constraints: 

for alH, k, = e ioo e oofc • (VD 

Now, when i 7^ 1 and k 7^ 1, the correlators that appear 
in (T7\ are already fixed according to (15[ , and the con- 
straint is indeed satisfied, as a consequence of the fact that 
for bilocal correlations, P(a,c\x,z) = P(a|a;)P(c|zj3 (or 
(A X C Z ) = (A X )(C Z )). When i = 1 and k ± 1, or vice 
versa, dTvT > gives linear constraints on the free correlators ejg^. 
and ejgj. There finally remains only one0 quadratic constraint 
in the case where i = k = 1: ejgj = ejoo eggi- 

Because the representation in terms of correlators nicely 
separates the fixed and free parameters of the (bi)local decom- 
position, it simplifies the search for explicit bilocal decompo- 
sitions quite significantly. We now present how this problem 
can be tackled in practice. 

3. Looking for explicit bilocal decompositions 

To determine whether a given correlation P is bilocal or 
not, one can now look whether there exist correlators e^j — 
some of which are fixed by (TT~5T > — that satisfy the nonnegativ- 
ity constraint ( TToT l and the bilocality assumption ( TTTb . 

This non-convex feasibility problem can be addressed from 
two different perspectives. First, one can try a heuristic search 
that will provide an explicit bilocal model if successful, thus 
proving P £ B — unless an exhaustive search can be under- 
gone, a negative result will however be inconclusive. Sec- 
ondly, one can try to solve a convex relaxation l23ll of the 
quadratic constraints (T7\ in the non-convex problem — a neg- 
ative result proving P ^ B. Both approaches lead naturally 
to the usage of numerical algorithms to solve the non-convex 
problem or its convex relaxation^. 

In practice, one may be interested in the robustness of a 
non-bilocal correlation to experimental imperfections, such as 
noise (see in particular subsections IIII A ll and llllBI ) or detec- 
tion inefficiencies (see subsection IIII C\ . One can then trans- 
form the above feasibility problem into an optimization prob- 
lem; for instance, one may want to maximize the noise that is 
tolerated by a correlation (or minimize the visibility of a noisy 
correlation of the form P(V) = VP+ (1 - V)P n . see eq. d33b 
below) before it becomes bilocal. The brute force search will 
give an upper bound on the tolerated noise (a lower bound on 



4 The fact that P(a,c\x,z) = P(a\x)P(c\z) directly follows from the bilo- 
cality condition 0: after summing over Bob's outputs 6, the two integrals 
over Ai and A2 factorize. Note that this equality also holds for quantum 
correlations established from independent sources, as in Fig.fTJ 

5 For more than 2 inputs for Alice and Charlie, there will remain more than 
one quadratic constraint, but the correlators representation will still signif- 
icantly simplify the search for explicit bilocal decompositions. 

6 In our numerical tests, we used the standard optimization toolbox of 
MATLAB, which implements the algorithms described in f24ll . and the 
BMIBNB solver of YALMIP for convex relaxations. 
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the visibility threshold VbUoc, as defined in subsection MI A II 
below), while the convex relaxation will give a lower bound 
on the tolerated noise (an upper bound on VbUoc)- 

Examples of explicit bilocal decompositions, for the quan- 
tum correlations studied in section [HI] are given in Ap- 
pendix [B] 



D. Non-linear inequalities for bilocal correlations 

As mentioned earlier, another approach to study the 
(non-)bilocality of given correlations is to test Bell-type in- 
equalities. Indeed, one can derive in some cases analytical 
constraints satisfied by all bilocal correlations; if a correlation 
is found to violate these constraints, then this implies that the 
correlation is non-bilocal. These non-linear bilocal inequali- 
ties thus provide convenient tests of non-bilocality, which are 
directly implementable in experimental demonstrations. 

We start by deriving such a bilocal inequality for the case 
where the three parties all have binary inputs and outputs. We 
then show that our first inequality implies similar inequali- 
ties for the case where Bob now has only one possible input, 
with 4 or 3 possible outputs (while Alice and Charlie still have 
binary inputs and outputs). These cases will be relevant for 
entanglement swapping experiments, as we will later see in 
section Hill 



1. A bilocal inequality for binary inputs and outputs 

We thus consider first the scenario where Alice, Bob and 
Charlie have binary inputs and outputs x, y, z, a, b, c <E {0, 1}. 

Let us define, for a given correlation P 22 , the tripartite cor- 
relation terms 

(A x B y C z ) P 22 = £ {-l) a+b+c P 22 {a,b,c\x,y,z) 

a,t>,c 

and the following linear combinations I 22 , J 22 : 

/22 = i £ (^BqOp^, (18) 

a;, 2=0,1 

J 22 = \ (-l) x+a (^SxC z )p 22 . (19) 

x,z=0,l 

As we show below, if P 22 is bilocal, then the following 
nonlinear inequality necessarily holds: 



Using the fact that |(-Bo)ai.a 2 I < 1, 



< - JJ dX 1 dX 2 p 1 (X 1 )p 2 (\ 2 )\(A ) Xl +(A 1 ) Xl \\(C )x 2 +(C 1 ) X2 

< dXipi(Xi) 1 '-x dX 2 p 2 (X 2 ) 1 1 



and one can show similarly, that 

|J 22 | 

< /dAxp^AijJ L x/dA2 / 92(A 2 ) J 



Now, for any r, s, r' , s' > 0, the inequality yjrs 
yjr + r'V 's + s' holds. Applied to the above two bounds on 
|/ 22 | and I J 22 |, we obtain 



< 



v 



dAipi(Ai) 



(A )x 1 +(A 1 )x 1 , (A )x-(A 1 )x 1 



dX 2 p 2 (X 2 ) | 



(Co)a 2 +(Ci)a 2 , (Co)a 2 — (Ci)a 2 



VW\ + \fW\ < 1 



(20) 



Furthermore, |(A ) Al + (Ai) Xl |/2 + \(A )x 1 - (Ai) Al |/2 - 
max(|(Ao)A 1 | ! |(Ai)x 1 |) < 1 and similarly, |(C )a 2 + 
(d) A2 |/2 + |(C )a 2 - (Ci) Aa |/2 < 1. After integrating 
over Ai and A2 in the previous expressions, we obtain inequal- 
ity ©. □ 

A projection of the correlation space onto the (J 22 , J 22 ) 
plane is shown on Figure [3] Note that the bilocal inequal- 
ity < f20b is tight in this plane, i.e., any values of I 22 and J 22 
such that \/\T^\ + \/\J 22 \ < 1 can be obtained by a bilocal 
correlation; see Table Q] in Appendix IB1 for an explicit bilocal 
decomposition. One can show on the other hand that the local 
correlations satisfy |/ 22 | + | J 22 | < 1 (recall that C is the con- 
vex hull of B), which can be understood as Bell inequalitiefl 
Finally, non-signaling correlations — such that their marginal 
probability distributions when one discards some parties do 
not depend on the settings of the discarded parties — also form 
apolytope, bounded by max(|/ 22 |JJ 22 |) < 1; in particular, a 
non-signaling correlation such thatj a + b + c = xy + yz (mod 
2), which can easily be realized with two Popescu-Rohrlich 
(PR) boxes 12611 (one shared by Alice and Bob, one shared by 
Bob and Charlie |H3l), reaches the values J 22 = J 22 = 1. 



Proof. By assumption, P has a bilocal decomposition of the 
form©. Defining {A x ) Xl = E a (-^) a P 22 (a\x, Ai), and 
with similar definitions for (B y )x 1 .x 2 an d (Cz)a 2 > one 8 ets 

I 22 = \[[d\ 1 d\ 2Pl (\ 1 )p 2 (\ 2 ) 

((A )x 1 +(A 1 )x 1 )(B )x 1 ,x 2 {(Co)x 2 +(C 1 )x 2 )- 



1 For instance, the inequality I 22 + J 22 < 1 can be written as (a quarter of) 
the sum of two equivalent facet inequalities of the tripartite local polytope 
(it is therefore not a facet itself): (A B Co + A1B1C1 + A\B Cq — 
A0.B1C1} < 2and{A 1 B C 1 +A B 1 Co + A B C 1 -A 1 B 1 C ) < 2, 
which are facets of the 'Class 3' type as defined in 

8 More precisely, P 22 (a, b, c\x, y, z) = |[1 + (-iy+b+c+xy+yzy 
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FIG. 3. Projection of the tripartite correlation space in the (J, J) 
plane, for I = J 22 or I 14 and J = J 22 or J 14 , as defined in COD- 
1 1 9b or J21U22b . for the "22-" and "14-"cases, resp. The non-convex 
bilocal set B is delimited by the four portions of parabolas, corre- 
sponding to the inequality \/|T|+ \/\J\ < 1 (eqs. J20b or d23t). It is 
included in the local set C delimited by the dashed lines, for which 
l-^l + \J\ _■ 1- The set Q of quantum correlations is also limited 
in this plane by |7| + \ J\ < 1. In the 22-case, the non-signaling 
polytope MS 22 is delimited by the outer dotted square, defined by 
max(|J 22 |, | J 22 1) < 1; in the 14-case, the projection of the non- 
signaling polytope A/iS 14 coincides with that of the local and quan- 
tum sets. 

The figure can also be understood as a 2-dimensional slice of the cor- 
relation space, containing the quantum correlation Pq = Pq 2 i34l 
or Pq 4 j30l , the fully random correlation Po = Pq 2 or Pq 4 , and the 
bilocal correlations Pi = Pf 2 or P} 4 and Pj = Pf or P} 4 , as 
defined in footnotes Q/3] and [T7] In this slice, the bilocal set is star- 
convex. One can see that the quantum correlation Pq = iPj + \Pj 
is local, but not bilocal. When adding some white noise, it enters the 
bilocal set for visibilities V < | (see subsection IIII A It . 



2. Scenario with one input, four outputs for Bob 

We now consider the case where Bob has only 1 possible 
input and 4 possible outputs, while Alice and Charlie still 
have binary inputs and outputs; this could in practice corre- 
spond to the case where Bob performs a complete Bell state 
measurement, in an entanglement swapping experiment — see 
section IIII Al below. 

Let us denote Bob's outputs by two bit^| b = bPb 1 = 
00,01,10 or 11, and by P 14 (a, c\x, z) the correlation 



shared by Alice, Bob and Charlie. Similarly to the previous 
case, we now define the tripartite correlation terms 

a,fc°bV 

and the linear combinations I 14 and J 14 as follows: 

I 14 = \ E (A*B°C Z } PI *, (21) 

it, 2=0,1 

J 14 = \ E {-1) X+ '{A«B 1 C,) P1 A. (22) 

X, 2=0,1 

As in the previous case, if P 14 is bilocal, then the following 
nonlinear inequality necessarily holds: 

VW\ + VW\ < 1- (23) 

Note that this implies in particular the inequality previously 
derived in lfl5ll (see appendix 0. 

Proof. We show that inequality (l23l l can directly be derived 
from (1201 1. Indeed, from the correlation P 14 , the three parties 
can obtain a correlation P 22 (a, b, c\x, y, z), with now binary 
inputs and outputs for Bob, if, for a given input y g {0, 1}, 
Bob simply outputs the corresponding bit b y . Formally: 

P 22 (a, &, c\x, y, z) = P 14 (a, b v = b, c\x, z) 

= J2 S b.byP 14 (a,b Q b\c\x 1 z). (24) 

One can easily check that for the correlation P 22 thus ob- 
tained, (A x B y C z )p22 = (A x B y C z )pn and the values of I 22 
and J 22 as defined in ( fT8ifT9l ) coincide with the values of I 14 
and J 14 obtained from <[2T}f22]). 

Suppose now that P 14 is bilocal. Since the processing from 
P 14 to P 22 is made locally by Bob, then P 22 is also bilocal, 
and therefore it satisfies (l20t . Since I 22 = I 14 and J 22 = 
J 14 , then (|23]) also holds. □ 

The projection of the correlation space onto the (I 14 , J 14 ) 
plane can also be seen on Figure [3] As before, the bilocal 
inequality d23l l is tight in this plane. On the other hand, local 
correlations satisfy the Bell inequalities |/ 14 | + | J 14 | < 1. 
Interestingly, it turns out that non-signaling correlations in this 
scenario are now also bounded b\Pl |/ 14 | + | J 14 | < 1. 

3. On the difference between the case with 1 input / 4 outputs 
and the case with 2 inputs / 2 outputs for Bob 

We observe a quite strong similarity between the cases 
where Bob has binary inputs and outputs (the "22-case"), and 



9 To clarify the notations, note that we use superscripts on Bob's output bits 
b° and 6 , to distinguish the case where they form one single output (b = 
b ^), from the previous case where bo and 6i were Bob's outputs for two 
different inputs (y = and 1, resp.). Note also that since Bob has only one 
possible input y, we do not need to specify it in P 14 (a, b^b 1 ,c\x, z). 



For instance, the non-signaling assumption implies I + J = 

sum is over all indices 

x, z, a, 6°6 1 , c such that 6° ^ a © c and b 1 ^ a © c © x ffi z, and 
where © denotes the addition modulo 2. 
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where he has 1 input with 4 possible outputs (the "14-case"). 
As we have just seen in the proof of d23l l, a scenario of the first 
kind can easily be obtained from a scenario of the second kind, 
and constraints on the correlations in the 22-case imply con- 
straints in the 14-case. However, this is only a one-way proce- 
dure. From a correlation with binary inputs and outputs, one 
cannot simply obtain a correlation with four outputs: indeed, 
bo and b\ in the binary case are incompatible measurement re- 
sults, and can in general not be outputted simultaneously to 
define a single 4- valued outcome b = bobi. 

From a more technical point of view, when looking for 
explicit (bi-)local decompositions, the definitions of weights 
q a py or of correlators ejjj for instance will be formally the 
same, the only difference being in the interpretation: in the 
22-case, /3q and (3i are different single-bit outputs, corre- 
sponding to different inputs, while in the 14-case, f3° and (3 1 
form a single 2-bit output, for a single input. Since the outputs 
(3° and /3 1 in the 14-case are compatible, the bilocal models 
are more constrained in that case than in the 22-case: correla- 
tors of the form ej n ^ (with i, k ^ 11) are fixed by correla- 
tions P 14 , but not by correlations P 22 . Although our bilocal 
inequality does not illustrate this fact, it might be the case that 
more severe constraints on bilocal correlations can be derived 
in the 14-case than in the 22-case. What can be seen how- 
ever when looking only at the linear combinations J 22 / 14 and 
j22/i4 j s tna ^ as a i rea( jy mentioned, the Bell inequalities of 
the form \I 22 \ + \J 22 \ < 1 for local correlations can be vi- 
olated by non-signaling correlations; but interestingly, their 
counterpart |P 4 | + | J 14 | < 1 cannot (see Figure[3]l. 



4. Scenario with one input, three outputs for Bob 

Let us finally consider the "13-case", where Bob has only 
1 input and 3 possible outputs (again, Alice and Charlie still 
have binary inputs and outputs); in practice, this could cor- 
respond to an incomplete Bell state measurement — see sec- 
tion |mB2]below. 

To compare with the case with four outputs for Bob, we 
still use two bits to denote Bob's outputs, b = bPb 1 = 00, 01 
or {10 or 11}: here, Bob does not distinguish his outcomes 
10 and 11, i.e. [tfb 1 = 10] = [b% x = 11]. We denote 
by P 13 (a, b, c\x, z) the correlation shared by Alice, Bob and 
Charlie. 

By analogy with the previous cases, let us now define the 
following tripartite correlators: 

(A x B a C z ) P i3 

= E (-1) Q+C [^ 13 (a, 00, c\x, z) + P 13 (a, 01, c\x, z) 

-P 13 (a, {10 or 11}, c|z,z)], 

(A X B C z ) pi3 !b o =0 

= E (- 1 ) a+c [ pl3 ( a > 00 ' c ^ z ) - pl3 ( a ' 01 » z )\ > 

a,c 

and, in a similar way again as before, the following linear 



combinations: 

/13 = i E (A x B a C z ) P is, (25) 

3!,Z=0,1 

J 13 = \ E (-ir+'^B'C^P" tf>=o ■ (26) 

27,2=0,1 

Once again, all bilocal correlations P 13 necessarily satisfy 



VWl + VWl < 1, (27) 
and this inequality is tight in the (I 13 , J 13 ) plane. 

Proof. We show here that inequality (|27T i can directly be de- 
rived from d23l . Indeed, from the correlation P 13 , the three 
parties can obtain a correlation P 14 (a, b, c\x, y, z), with now 
four possible outputs for Bob, in the following very simple 
way: when Bob gets an outcome b = 0b 1 , he outputs it di- 
rectly; when he gets the outcome b = {10 or 11}, he outputs 
b = 10 or b = 11 at random. Formally: 

P 14 (a,06\ckz) = P 13 (a,0b\c\x,z) 

(28) 

P u (a,lb 1 ,c\x,z) = iP 13 (a, {10 or 11}, c\x,z) . 

One can again check that for the correlation P 14 thus ob- 
tained, (A x B°C z ) P i4 = (A x B°C z ) P i3 and(A x P 1 C 2 ) P i4 = 
(A X B 1 C Z ) pi3(,o_ ; therefore the values of J 14 and J 14 as de- 
fined in ( T2T1 - |22] | coincide with the values of J 13 and J 13 ob- 
tained from d25l - [26] >. 

Suppose now that P 13 is bilocal. Then so is P 14 , which 
therefore satisfies ((23j. Since I 14 = I 13 and J 14 = J 13 , then 
P 13 satisfies (27). □ 

The projection of the correlation space onto the (I 13 , J 13 ) 
plane is shown on Figure [4] As in the 14-case, the lo- 
cal and non-signaling sets are delimited in this plane by 
|/ 13 | + | J 13 1 < 1, while the bilocal correlations satisfy 
-\/|/ 13 | + y/\ J 13 1 < 1 (which is tight in this plane: see Ta- 
bleHUin AppendixlBj. 

It is also relevant to restrict ourselves to correlations with 
a random marginals (P°)pi3 = for Bob (with an obvious 
notation) and random bipartite marginals (A X C Z ) pi3 = for 
Alice-Charlie, such as the quantum correlation Pq 3 ( T3"6i l stud- 
ied in subsection HUB 21 below. With these additional con- 
straints, the local and non-signaling sets are delimited in the 
(J 13 , J 13 ) plane b£] |P 3 | + 2| J 13 | < 1; see Figure g] 



The non-signaling assumption indeed implies for instance I 13 + 2 J 13 = 
1 + (B°) -jJ2x z( A * C z) - 2 E pl3 ( a .°' ,1 . c l a: :- 2 ). where the last 
sum is over all indices x, z, a, b 1 , c such that a ^ c and b 1 = x ® z. If 
Bob's marginal (B°) and all Alice-Charlie's bipartite marginals (A X C Z ) 
are zero, it follows that i" 13 + 2 J 13 < 1. 



9 




-1 -1/2 1/2 1 



FIG. 4. Projection of the tripartite correlation space in the (J 1 , J ) 
plane, as defined in | |25H26| |. for the "13-case". The projections of 
the bilocal set B 13 , of the local polytope £ 13 , of the quantum set 
Q 13 and of the non-signaling polytope A/5 are similar to those of 
Figure [3] When restricting to correlations with a random marginal 
(B°} P i3 for Bob and random bipartite marginals (A X C Z ) pis for 
Alice-Charlie, the local (£rn D ), quantum (Qrnd) ar, d non-signaling 
(A/S^nd) sets are instead delimited by |7 13 | + 2| J 13 < 1 (dashed 
diamond). 

The restriction of the figure to the dashed diamond can also be un- 
derstood as a 2-dimensional slice of the correlation space, containing 
the quantum correlation Pq 3 i36l , the random correlation Pq 3 , and 
the bilocal correlations P} 3 and P} 3 , as defined in footnote[T8] One 
can see that the quantum correlation Pq 3 = §ij 3 + ^Pj 3 is local, 
but not bilocal. When adding some noise, it enters the bilocal set for 
visibilities V < | (see subsection IIII B 21 . 



III. QUANTUM VIOLATIONS OF BILOCALITY IN 
ENTANGLEMENT SWAPPING EXPERIMENTS 

Since John Bell's work in the 1960's fl], it is well under- 
stood that his locality assumption can be falsified by quan- 
tum correlations. Hence, our bilocality assumption can a for- 
tiori also be falsified quantum mechanically. Since the latter 
is a stronger assumption than the former, one may wonder 
whether it can lead to stronger tests of quantumness; we show 
now that this is indeed the case. 

To justify the use of the bilocality assumption, we consider 
below scenarios where two independent quantum sources Si 
and S*2 send particles to Alice and Bob in the state Qi, and 
to Bob and Charlie in the state Qi, respectively, so that the 
overall quantum state is 

QABC = Ql ® 02 ■ (29) 

This typically corresponds to entanglement swapping experi- 
ments ||l 1J, as depicted in Figure Q] 

We always assume below that Alice and Charlie have bi- 
nary inputs and outputs. As for Bob, we consider the case 
where he can perform a full Bell state measurement (he has 
only one possible input, with four possible outputs), the case 
where his measurement results group the Bell states two by 
two, and he can choose among two possible pairings (he has 
binary inputs and outputs), and the case where his partial Bell 
state measurement distinghishes two Bell states, but groups 
the other two together (he has again only one possible input, 
with now three possible outputs). 

In most cases (with the notable exception of subsec- 
tion IIIIC 21 below), our bilocal inequalities derived in the 
previous section will be sufficient to demonstrate the non- 
bilocality of the quantum correlations thus obtained. 



5. Do our bilocal inequalities fully characterize the bilocal set? 

Note that one can of course also derive many equivalent 
versions of inequalities (|20| |. ( l23l l and (l27l i, where the in- 
puts and/or outputs are permuted. One may wonder whether 
these inequalities are enough to delimit the bilocal set, as it 
is the case for instance with the Clauser-Horne-Shimony-Holt 
(CHSH) inequality lElll and its equivalent versions for 2 par- 
ties with binary inputs and outputs, which — together with triv- 
ial inequalities of the form P > — fully characterize the cor- 
responding local set 1I22I1 . The answer is negative: there exist 
nonbilocal correlations that satisfy (l23l l for instance and all 
its equivalent versions. Although in most practical cases we 
study in the next sections, considering inequalities (l20l l. (l23l l 
or (ITTl i will be enough to demonstrate non-bilocality, the study 
of the detection loophole in section [III C I will provide an ex- 
ample where these are not sufficient; we will then resort to 
convex relaxation methods. 

Fully characterizing the bilocal set, in the scenarios con- 
sidered here, by a simple list of inequalities remains an open 
problem. 



A. Entanglement swapping experiment 
with a complete Bell state measurement 

We start by considering a typical standard entanglement 
swapping experiment, where the sources Si and S2 each pro- 
duces a singlet state l^ - ), and where Bob performs a com- 
plete Bell state measurement on the two particles he receives 
from the two sources; the four possible outcomes b = bPb 1 = 
00,01,10 or 11 he can obtain correspond to the four Bell 
states (with standard notations) |<I> + ), |$~), \^> + ) or |* _ ), 
respectively. Alice and Charlie can each choose a projec- 
tive measurement (with binary outcomes) to perform on their 
qubit, described by the observables A x and C z (corresponding 
to their inputs x and z, resp.). 

Suppose that Alice and Charlie can measure either Aq = 
C = (<7 Z + v x )/V2 (for x, z = 0) or i x = C x = (a z - 
fix)/ V2 (for x, z = 1) on their particle, where a z and <j k are 
the Pauli matrices^- Quantum Mechanics predicts that Alice, 



In order not to confuse the notations, we use the fonts X, Y, Z for the Pauli 
matrices <r x , ov and ct z or for directions on the Bloch sphere, and the fonts 
x, y, z for Alice, Bob and Charlie's inputs. 
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Bob and Charlie will observe the following correlation: 

(-l) b ° + (— ■ 



P Q *(a,b%\c\x,z) 
1 

" 16 



1 + (-1) Q+C 



From the definitions d2Tl - l22l i. one easily obtains 

W) = J"(? Q ") = i, 



(30) 



(31) 



which violates d23l l. One can thus conclude that the quantum 
mechanical correlation Pq 4 (f3Qb is non-bilocal. 

It turns out however that Pq 4 is locaO Like local and non- 
signaling correlations, quantum correlations are bound to sat- 
isf>0 |/ 14 | + | J 14 1 < 1 (see Figure EJ; note that thisjiolds 
even if the state qabc does not have the product form I 



1. Resistance to noise 

An interesting figure of merit to quantify the 
non-(bi-)locality of a correlation is its resistance to noise. 
One way to model noise£3 is to suppose that each source Si 
introduces white noise with probability 1 — Vi, i.e. corre- 
sponding to a visibility instead of sending a pure (say, 
2-qubit) state \ipi), the state it actually sends is 



Qi(Vi)=V i \lp i )(^i\ + (l-Vi)l/4. 



(32) 



In the case we consider here, with maximally entangled 
states and random marginal probability distributions, the re- 
sulting quantum correlation will only depend on the product 



13 ^Q 4 can i n d ee d be decomposed as Pq 4 = ^Pj 4 + \P\ 4 ( see 

Figure[3j, with P} 4 (a, b ^ , c\x, z) = i[l + and 

P} 4 (a,&°& 1 1 c|x,.s) = ±{1 + {-l) x + z+a + c+b1 ]. P} 4 and P} 4 
are bilocal (and hence, local): they can be obtained from the explicit 
decompositions of Table [II] in Appendix [B] for (7 14 =1, J 14 =0) and 
(I 14 =0, J 14 =l), respectively. 

Note that Pj 4 and Pj 4 , as well as the fully random correlation 
Pq 4 , are all invariant with respect to the symmetries (a, b ,^) o 
(oel,6°el,6 1 ©l), (6°,6 1 ,c) «-> (6°el,6 1 ©l,c©l), (z.fc 1 ) o 
(zffil, ft 1 ©!) and (z, b 1 ) -H- (zffil, fc 1 ©!), which can all be applied bilo- 
cally (i.e., independently between A—B and B-C). When applying each 
of these symmetries with probability 1, any correlation P 14 (giving values 
I 14 , J 14 ) is projected onto a correlation P] 4 = I 14 P} 4 + J 14 Pj 4 + 
(1— J 14 — J 14 )Pq 14 on the 2-dimensional slice represented on Figure [5] 
This "depolarisation" is similar to that introduced in Ref. flUl ; similar de- 
polarisation processes can be defined in the other (22- and 13-) cases. 

14 Note that all values of I 14 , J 14 such that |/ 14 | + | J 14 | < 1 can be ob- 
tained quantum mechanically; for instance, for Alice and Charlie's mea- 
surement settings of the form Aq/i =C*o/i = cos ^zi s i n "ox , we get 



cos 

(with < I 14 , J 14 < 1) is recovered 



sin 9, and the full line segment I + J 



1 



15 Note that one could also consider some additional white noise due to the 
measurement apparatuses, by introducing a visibility tip for each party 
P = A, B or C. In the case we consider here, the resulting quantum 
correlation would then depend on the product V = v \v^v j\Vb v C °f the 
visibilities of each source and measurement apparatus, and our argument 
would remain unchanged. 



V = V1V2 of the visibilities of each source, and is simply 
given by 



14 



(33) 



where P 14 is the fully random probability distribution (i.e., 
P^ 4 {a,b°b\c\x,z) = 1/16 for all a, c, x, z). The 
largest visibility for which the correlation Pq 4 (V) admits a 
bilocal decomposition defines the bilocal visibility threshold 
Vbiioc, and can be used to quantify the non-bilocality of Pq 4 : 
the smallest Vuioc is, the more resistant to noise, and hence 
the more bilocal Pq 4 is. 

Noting that I 14 [Pq(V)} = J 14 [Pq(V)} = \V, we con- 
clude from our bilocal inequality (|23l l that Pq 4 (V) is non- 
bilocal for all visibilities V > 50% (see Figure [3]). On the 
other hand, for visibilities V < 50%, one can find an ex- 
plicit decomposition that proves that P^fiV) is bilocal (see 
Table [II] in Appendix [B); our inequality thus detects opti- 
mally the resistance to noise of the correlation Pq 4 , for which 
V bllo c = 50%. 

The visibility V can be understood as the visibility of the 
maximally entangled state that results from the entanglement 
swapping process. In order to check the nonlocality of that 
state in the standard locality scenario, one could test the 
CHSH inequality 12811 : this would require the use of different 
measurement settings for Alice or Charlie, and would require 
a visibility V > ^= ~ 70.7% for the CHSH inequality to 
be violated. Actually, no Bell inequality can be violated (us- 
ing Von Neumann measurements) for visibilities smaller than 
V ~ 66% l29ll . Our assumption allows one however to exhibit 
non-bilocal correlations for visibilities as low as 50%. This 
illustrates the advantage of the bilocality assumption, which 
simplifies the requirements for the demonstration of quantum- 
ness in entanglement swapping experiments l30l[3lll . 



2. On our choice of measurement settings 

The measurement settings we chose for Alice and Charlie 
above (and which were already introduced before in lfl5ll ) are 
the ones giving the bilocal quantum correlation Pq 4 (f30b that 
is the most resistant to noise we could find, i.e. with the small- 
est bilocal visibility threshold Vbiioc- 

These measurement settings were first obtained numeri- 
cally, following the approach introduced in subsection III CI 
Our extensive numerical tests convince us that we have found 
the optimal settings for the case where the sources send two 
singlet states, where Alice and Charlie can choose among two 
projective measurements, and where Bob performs a com- 
plete Bell state measurement. The symmetries of the quan- 
tum correlation Pq 4 < f30b then actually inspired our defini- 
tions ( F2T1 - I221 ) for P 4 and J 14 , and the whole analysis of sub- 
section HID I interestingly, our bilocal inequality d23l ) is suf- 
ficient to demonstrate the non-bilocality of Pq 4 (V) down to 

V = Vbiioc- 
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B. Partial Bell state measurements 

An ideal entanglement swapping experiment requires Bob 
to perform a complete Bell state measurement. This might not 
be a trivial thing to do; actually, it is known to be impossible 
to perform this ideal joint measurement with linear quantum 

optics m. 

From an experimental perspective, it is therefore interest- 
ing to study the consequences of the bilocality assumption in 
scenarios where Bob does not perform a complete Bell state 
measurement, but only a partial one. We thus consider be- 
low cases where Bob's measurement does not allow him to 
discriminate the four Bell states, but only subsets of the Bell 
states. 

For instance, he may perform a measurement that allows 
him to discriminate one Bell states vs the other three. In that 
case however, we found no advantage with the bilocality as- 
sumption, over a test of standard locality (for instance, a test 
of CHSH between Alice and Charlie, conditioned on Bob hav- 
ing observed the Bell state he can distinguish), even if Bob 
may have different possible inputs that allow him to choose 
which of the four Bell states he wants to discriminate: we 
always found V/ oc = V&«oc = V V%- Another possibility 
would be for Bob to distinguish pairs of Bell states. If he can 
only distinguish two states vs the other two, the correlation 
shared by the three parties will be bilocaQ if ne can choose 
among two different pairwise groupings of the same four Bell 
states (when he has two possible inputs), they will obtain non- 
bilocal correlations, with a visibility threshold Vbiioc = 50% 
(see subsection 1MB II below): the case where he can choose 
among the 3 possible pairwise groupings does not provide any 
advantage over the case with two inputs. Finally, Bob may be 
able to discriminate two Bell states perfectly, but not to dis- 
tinguish the other two. In that case, the correlation will again 
be non-bilocal, with now a visibility threshold Vuioc = 2/3 
(subsection IIIIB 2l below); when Bob has more than one in- 
put that allow him to choose which two Bell states he wants to 
discriminate perfectly, the situation is the same as the previous 
one, with pairwise groupings. 

1. Binary inputs and outputs for Bob 

Let us thus start by considering the case where Bob has two 
possible inputs, and he wants to distinguish either the l^*) 
vs the Bell states, or the |$ + /*+) vs the | <I> / ^ > 
Bell states. That is, he measures either Bq = |<f> + )($ + | + 
|$-)($-| - |# + )(#+| - |#-)(*-| = <7 Z ® o-z or Bi = 
|*+)(*+|-|$-)(*-| + |tf + )(tf+|-|*-)<tf-| = <7 X ®(T X ; 
note that his measurement is actually a separable measure- 
ment. 

We still assume that the two sources produce singlet states, 
and that Alice and Charlie perform the same measurements as 



This is due to the fact that for a pairwise grouping of the Bell states, Bob's 
measurement is separable, of the form try ® <ttj, with U = X, Y, or Z 
depending on the grouping (see llll BTl . 



before, in subsection llll Al Arj/i = Qi/i = (5z ± 6x)/ v2. In 
this scenario where all parties have binary inputs and outputs, 
the quantum correlation shared by the three parties is then 

P 22 (a, b, c\x, y, z) = I [1 + l(_ ir+ b +c+xy+y ^ ) (34) 
for which, from the definitions ( fT8ilT9l l, one gets 

7 22 (P 22 ) = J 22 (P 22 ) = \- 05) 

This violates d20l ), which proves that Pq 2 is not bilocal. 

The bilocal visibility threshold can be defined as in the 
previous section, by considering correlations of the same 
form as d33l (with now Prj 2 (a, b, c\x, y, z) = 1/8 for all 
a, 6, c, x, y, z). We find that Pq 2 (V) violates inequality d20b 
for visibilities V > 50%, and that it admits a bilocal model 
for visibilities V (see Table [I] in Appendix IB1. so that again, 
for P 22 , Vuioc = 50% (see Figure^. 

The practical consequences of this result are however not as 
interesting as in the 14-case. Bob's measurement being sepa- 
rable, we are not considering here an entanglement swapping 
experiment. In fact, the scenario here amounts to perform- 
ing two tests of the CHSH inequality, between Alice-Bob and 
Bob-Charlie. The requirement V > 50% simply corresponds 
to the requirement that at least one of the visibilities Vi of the 
CHSH tests be larger than 1 / \/2; however, this is already a 
sufficient condition to demonstrate simply Bell nonlocality. 

For completeness, note that as was PA 4 in the 14-case, Pq 2 
is locaQ 

Like local correlations (but unlike non-signaling 
ones, here), and as in the 14-case once again, quantum corre- 
lations actually also satisfy |/ 22 | + | J 22 | < 1, which is again 
tight and holds even if the state qabc does not have the prod- 
uct form d29l ); this can be seen by expanding the factors of the 
positive operator 6 = (A + - B C + ) 2 + (C + - A+B Q ) 2 + 
(i_ - Pi(7_) 2 + (6- - i_Pi) 2 , with i± = A ± ii and 
C± = Cq±Ci (and where tensor products are implicit, while 
identity operators are omitted; e.g., BqC+ actually stands for 
1 g> B ® C+). 

2. One input, three possible outputs for Bob 

Another interesting case is when Bob performs a single in- 
complete Bell state measurement, with now three outcomes. 
We consider here for instance the case where his outcomes 
b = b^b 1 = 00, 01, {10 or 11} correspond, respectively, to 
| ( E )+ ), l^ - ) and l^*) (whichBob cannot discriminate, as they 
give the same outcome). Such a measurement can be realized 



P'q can indeed be decomposed as Pq 2 = \Pf 1 + ^Pj 2 ( see 
Figure fj), with Pf 2 (a, b, c\x, y, z) = ±[1 + 5y,o(-l) a+b+c ] and 
P] 2 (a,b,c\x,y,z) = ^[1 + S yA {-l) x + z + a+b + c }. Pf 2 and Pf are 
bilocal: they can be obtained from the explicit decompositions of TableU 
in Appendix|B] for (i" 22 = l, J 22 =0) and (7 22 =0, J 22 =l), respectively. 
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with linear quantum optics rf32"tl . hence the practical motiva- 
tion for studying this particular case. 

We consider now the following measurements for Alice and 
Charlie: A = C = (V2a z + a x )/V3 and A x = C x = 
(y/2a z — <7x)/ Vj. The quantum correlation Pg 3 (a, b, c\x, z) 
shared by Alice, Bob and Charlie is then 

pgia, ob\ c \ x , z) = m*+ (-ir +c 2+{ -T^ ] (36) 

Pg(a, {10 or 11}, c\x, z) = \ [l - |(-l) a+c ] 
One easily obtains, from the definitions (l25l - l26T i. 

r 13 (Ph 3 ) = l> J 13 (Pq 3 ) = 1, (37) 

which violates ( 127) . thus proving that Pq 3 is not bilocal. 

For noisy correlations of the same form as Q31 l. with now 
P 13 (a,06 1 ,c|a;,z) = ± and P 13 (a, {10 or 11}, c\x,z) = 
|, the above values of J 13 and J 13 are again simply to be 
multiplied by V. They violate inequality (|27| | for all V > §. 
For V < | on the other hand, one can find an explicit bilocal 
decomposition for Pq 3 (V), as in Table iHll of Appendix iBl 
Hence, the bilocal visibility threshold of Pq 3 is Vuioc = § 
(see Figure|4]i. 

Even when a complete Bell state measurement is not pos- 
sible, the bilocality assumption thus provides an advantage — 
compared to the standard Bell locality assumption — for prac- 
tical demonstrations of quamtumness in entanglement swap- 
ping experiments. 

Let us finally note that as in the previous cases, Pq 3 is lo- 
Like local and non-signaling correlations, quantum cor- 
relations with random single- and bi-partite marginals (as ob- 
tained from singlet states and a partial Bell state measurement 
for Bob) satisf)B |/ 13 | + 2\J 13 \ < 1; see Figureg] 

C. Resistance to detection inefficiencies 

Another experimental imperfection that is important to take 
into account is the fact that Alice, Bob and Charlie's detec- 
tors might not be 100% efficient. In a typical demonstration 
of nonlocality, this may open the well-known detection loop- 
hole 13311 . if the parties post-select their correlations on de- 
tected events only. 



Pq 3 can indeed be decomposed as Pq 3 = §P/ 3 + ^Pj 3 
(see Figure D, with P} 3 {a, Ob] c\x, z) = i[l + (-l)°+ c ], 
P} 3 (a,{10 or 11}, c\x, z) = |[l-(-l) a+c ] and Pj 3 (a, Ob 1 , c\x, z) = 
^[1 + (_l)*+z+a+c+'> 1 ], P 13 (a, {10 or 11}, c\x,z) = §. P} 3 and 
Pj 3 are bilocal: they can be obtained from the explicit decompositions of 
Table fTTTl in Appendix E for (J 13 =l, J 13 =0, K 13 =l, L 13 =M 13 =0) 
and (I 13 =0,J 13 = \,K 13 = -I, L 13 =M 13 =0), respectively. 
Note again that all values oil 13 , J 13 such that \I 13 1 + 2| J 13 1 < 1 can be 
obtained quantum mechanically; for instance, for the same measurement 
settings as in footnote [14] the full line segment I 13 + 2 J 13 = 1 (with 
< / 13 ,2J 13 < 1) is recovered. 



Restricted classes of local models with independent sources 
were actually considered before in lfl3l [HI] , and were pre- 
cisely studied in the context of the detection loophole. Here 
we initiate the study of the detection loophole with respect to 
the general assumption of bilocality, by considering the sim- 
plest cases, where only one party has an limited detection ef- 
ficiency r\ e [0, 1] while other parties have perfect detectors, 
and the case where both Alice and Charlie have the same de- 
tection efficiency r\ while Bob has 100% efficient detectors. 
More complex cases, for instance where all parties may have 
non-prefect detectors, are beyond the scope of this paper, and 
are left for future work ll34ll . 

In the preliminary study below, we consider again binary 
inputs and outputs for Alice and Charlie, and a complete Bell 
state measurement for Bob. We assume that when one party 
fails to get a conclusive result, they still output a result from 
their standard set of possible outcomes, either at random or 
according to a specific strategy — note that when the parties 
output random results, the situation is the same as for white 
noise in their measurement apparatus, and as that of imperfect 
visibilities studied before. The case where a no-detection re- 
sult is explicitly treated as a different outcome is left for future 
work; it is an open question, whether this may increase here 
the resistance to detection inefficiencies. 



1. Only one party has inefficient detectors 

We start with the case where one party (either Alice, Bob 
or Charlie) has imperfect detectors with efficiency rj, while 
the other two have perfectly efficient detectors. Note that in 
the case of Bob's Bell state measurement, if one of his de- 
tectors does not click, he might still get some partial infor- 
mation (such as on a subset of possible Bell states); we as- 
sume however that he does not make use of that information — 
which could possibly lead to a better resistance to detection 
inefficiencies — and we leave this potential improvement as an 
open research problem. 

We found that the best strategy is for Alice and Charlie 
to use the same measurement settings as in subsection 1111 Al 
which give the correlation Pq 1 d30l ). and for the party with 
inefficient detectors to output a random result in case of a no- 
detection event; in that case, similar to that studied previously, 
the resulting correlation is non-bilocal for rj > ribiioc = 50% 
(where the bilocal detection efficiency threshold r\uioc is de- 
fined in a similar way as the bilocal visibility threshold Vbiioc 
above); our bilocal inequality ( 1231 detects optimally its non- 
bilocality. 



2. Alice and Charlie both have imperfect detectors 
(Va = rjc = T),r]B = 1) 

The second case we consider is that where Alice and Char- 
lie both have imperfect detectors, with the same efficiency 
rj, while Bob has perfect detectors. Alice and Charlie still 
perform the same measurements as before. In case of non- 
detections, the best strategy is for Alice to output her input 



13 



67% 71% 75% 



Bilocal correlations 
Upper bound from our 
bilocal inequality 



67% 71% 75% 

Detection efficiency r] 



FIG. 5. (Non-)bilocality of correlations obtained by the strategy de- 
scribed in the main text, in the case where Alice and Charlie's de- 
tectors have a detection efficiency 77. The correlations are bilocal 
for visibilities V = V1V2 < V b v iloc , w ' m ^blioc depending on 77. 
As explained in the text, V biloc was estimated by matching lower 
bounds given by explicit bilocal decompositions (see Table llVl of Ap- 
pendix|B) and upper bounds obtained from convex relaxations of the 
bilocality constraint; our inequality j23\ also gives an upper bound 



on V b v iloc ; this bound was found to be tight only for 77 > 



directly (i.e., a = x), and for Charlie to always output c = 0. 
The resulting correlation is found to be bilocal for all 77 < 2/3; 
an explicit bilocal decomposition is given in TablellVlof Ap- 
pendix [B] On the other hand, convex relaxation techniques 
mentioned in subsection III C 31 allowed us to establish a nu- 
merical upper bound of 77 = 2/3 + e, with e ps 10~ 6 , above 
which the correlation is non-bilocal. We conclude that in this 
case, rjbiioc = 2/3. 

Taking also into account the noise in the state preparation, 
as in the previous subsections, the correlation then depends on 
77 and on V = v\V2- For any fixed value of 77, one can estimate 
the corresponding bilocal visibility threshold V ba , as shown 
on Figure [5] This was obtained again by comparing a lower 
bound on V h v iloc given by the explicit bilocal decomposition 
of Table [TV] in Appendix IB"l (to which we refer for analytical 
expressions for Vy, ), with a numerical upper bound derived 
using convex relaxations; the two bounds match again up to 
e ps 1CT 6 . 



Note that our bilocal inequality ( l23l can also be used to 
obtain an upper bound on ?7m oc and V biloc . However, this 
upper bound on VZL is found to be tight only for 77 > 3/4 
(see Figure [5]); and for V = 1, inequality d23l is violated 
only for 77 > 1/V2. This illustrates the fact that our bilocal 
inequality is not always sufficient to detect the non-bilocality 
of a correlation (see subsection lHD 51 ). 



3. Open problems related to the detection loophole 

A more complete study would be necessary to draw any 
conclusion regarding the advantage of the bilocality assump- 
tion compared to the locality assumption, with respect with 
the detection loophole. 

One should in particular consider the case where all three 
parties have inefficient detectors; note that because their mea- 
surements are inherently different, there is no reason to as- 
sume that Alice and Charlie should have the same detection 
efficiencies as Bob; rather, it would be relevant to consider 
a practical (and incomplete) Bell state measurement and see 
how Bob's efficiency would compare to that of Alice and 
Charlie Q. 

Many questions are left open here, including whether using 
partially entangled or higher-dimensional states may lower the 
required efficiencies, as it is the case in the standard scenario 
of Bell nonlocality SHI]. 



IV. FURTHER ISSUES ON QUANTUM NON-BILOCALITY 

We now present some additional results related to the study 
of quantum non-bilocality, coming back to the case where Bob 
performs a complete Bell state measurement and where all 
parties have perfect detectors. We first study the relation be- 
tween resistance to noise with respect to locality and to bilo- 
cality. We then present quantum violations of bilocality using 
partially entangled qubit states. Finally, we address the ques- 
tion of simulating quantum correlations in a bilocal manner. 



A. Trade-off between Bell nonlocality and non-bilocality 
for quantum correlations 

We have shown that the quantum correlations Pq 4 ( l30t , 
Pq 2 (f34t and Pq 3 d36t we exhibited are not bilocal. They are 
in fact the most robust to noise we could find in each scenario, 
and become bilocal for visibilities smaller than Vuioc = 50% 
in the first two cases, and V b u c = 2/3 in the last case. 

However, these correlations were found to be local. In each 
case, one can also, of course, obtain nonlocal quantum corre- 
lations by rotating for instance the measurement settings of 
Alice and/or Charlie. The nonlocality of such correlations 
can also be quantified by their resistance to noise: i.e., for a 
given correlation one can define the local visibility threshold 
Vioc below which the corresponding noisy correlation of the 
form (l33l becomes local. One may then wonder how the two 
visibility thresholds, Vi oc and V b u oc , behave, one compared to 
the other, when the measurement settings of Alice and Charlie 
vary. 

To illustrate the trade-off between the local and bilocal vis- 
ibility thresholds, let us consider here the first scenario (the 
14-case), where Bob performs only one measurement, with 
four possible outcomes. For the correlation d30l ), we had 
Vioc = 1 and VbUoc = \- By changing the measurement set- 
tings Aq,Ai,Cq and C\, while still considering a complete 
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trade-off between optimal resistance to noise 
of nonlocality and of non-bilocality 



Local visibility threshold Vi oc 

FIG. 6. Shaded area: local (V; oc ) versus bilocal (Vuioc) visibility 
thresholds for quantum correlations obtainable in standard entangle- 
ment swapping experiments, with binary inputs and outputs for Alice 
and Charlie. For the quantum correlation Pq 4 of eq. J30b for instance, 
on has (V ioc , Vuioc) = (1, |)- 



Bell state measurement for Bob, we modify the quantum cor- 
relation Pq 4 and thus obtain different corresponding pairs of 
visibility thresholds (Vi oc , Vbiioc)- From a numerical investi- 
gation, we obtained the set of pairs (Vi oc , Vbu oc ) represented 
on Figure [6] 

The upper boundary illustrates the fact that, obviously, 
Vbiioc < Vioc, note that the extreme left point for which 
Vioc = Vbiioc = -7= — 70.7% can be obtained for instance 
when one uses the standard settings to test the CHSH inequal- 
ity j28ll between Alice and Charlie (e.g., <r z and a x for Alice, 
(<7 Z ± Ox)/ a/2 for Charlie). One can see however that Vbiioc 
can be lowered down to Vbiioc — 2 — V2 ~ 58.6% while still 
having Vi oc = -j=. 

The lower boundary is of particular interest, as it expresses 
the trade-off between optimal resistance to noise of non- 
locality and of non-bilocality; it is the result of a multi- 
objective optimization problem. This front can for instance 
be parametrized by considering the following measurement 
settings: 



A = Co = cos 0q rr z + sin 6$ <7 X , 
A\ = C\ = cos Of fj z + sin Of a x , 

with 0f = (-l)*f-£f, and £e [0,1]. 



(38) 



A straightforward calculation gives, for the quantum corre- 
lation -Pq 4 (£) thus obtained and for the definitions (l2"Dj2~2"l >. 
I 14 = J 14 = |[1 + cos(£?)]. For a given visibility V, one 
just has to multiply these values of J 14 and J 14 by V; the cor- 
relation then violates inequality d23l l for all visibilities V > 
1/[1 + cos(£f )]. For visibilities lower than 1/[1 + cos(£f )] 
on the other hand, one can find an explicit bilocal decompo- 
sition, given in Table [V] of Appendix iBl Hence, this value is 



precisely the bilocal visibility threshold V biloc for the quantum 
correlation obtained with the measurement settings ( 1381 ). 

To test for locality, we note that for a tripartite correlation 
to be local, a necessary condition is that the corresponding 
bipartite correlation between Alice and Charlie, conditioned 
on one particular result of Bob, is local l37h . and therefore (in 
our case where Alice and Bob have binary inputs and outputs) 
it must satisfy the CHSH inequality 12811 . When conditioned 
on Bob obtaining \&~), we find that the value of the CHSH 
polynomial is CHSH = 2[cos(£ : |) + sin(£-j)]. For a given 
visibility V, the value of CHSH is simply multiplied by V, 
and the CHSH inequality CHSH < 2 is violated for all V > 
l/fcos^f) + sin(£ j)]. On the other hand, one can check 
that the corresponding correlation is local for visibilities lower 
than l/[cos(c; : |)+sin(c; f)]; an explicit local decomposition is 
given in TablefVlof AppendixlBl Hence, this value is precisely 
the local visibility threshold Vi oc for the quantum correlation 
we consider. 

The lower front in Figure [6] can thus be parametrized as 



V, 



bill. 



cos(£f) + sin(£f)'l + cos(£f) 



(39) 



with £ e [0, 1]. For £ = in particular, one gets the most 
non-bilocal correlations we could find, namely Pq 4 as in < T30b ; 

for £ = 1, we obtain the point (V^f , V b \= \) = 2 - V2) 
that we briefly mentioned before. 

We also looked at this relation between the locality and 
bilocality visibility thresholds in other scenarios and found 
similar trade-offs, although with different quantitative results. 
For instance, in the 22-case we found numerically that when 
Alice and Charlie's settings vary, while Bob's measurements 
are fixed, all correlations for which Vi oc = 1/ v2 also have 
Vuioc = l/v2 = Vioc', an improvement in Vuioc can only be 
obtained by increasing VJ oc . 



B. Non-maximally entangled states 

One might wonder whether our inequality d23l for bilocal- 
ity can be violated by non-maximally entangled states, and 
how non-bilocal the resulting correlations can be. 

Let us thus consider the case where the sources Si and S2 
send 2-qubit entangled states of the form 



cos^lOl) 



sin -1 1 10) 



(40) 



with 9i E [0, -|] (and i = 1,2). We assume that Bob performs 
a complete Bell state measurement (in the standard Bell ba- 
sis), and that Alice and Charlie can each choose among two 
projective measurements to perform on their respective qubits. 
The tripartite correlators that appear in the definitions (f2TT - l22b 
of I 14 and J 14 are easily found to be 



(A x B a C z ) P i, 
(A X B 1 C Z ) pi4 



a^c^ sin 6 1 sin 9 2 



(41) 
(42) 



where a| x and cf x are the Z and X components of the vectors 
a x and c z representing Alice and Charlie's measurements in 
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the Bloch sphere, respectively (for inputs x and z). We thus 



obtain 



i 14 = \K 
J 14 = 4(^ 



■ c?) sin 9i sin 6*2 . 



(43) 
(44) 



One can easily see that, in order to maximize y/\I\ + y \J\, 
the optimal settings of Alice and Charlie should be in the ZX 
plane, symmetric around the Z axis. More precisely, we find 
that the optimal measurements are, for both Alice and Charlie, 



<j z ± V sin 0i sin 62 <7 X 
V^l + sin^! sin 6 2 



leading to 



rl4 



1 



1 + sin 6*i sin( 



J 
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sin 2 6\ sin 2 62 
1 + sin 6\ sin 6*2 ' 



and 



sm f 1 sin #2 



(45) 



(46) 



(47) 



If sin6»i sin6» 2 > 0, we have ^f\I u \ + \/\J 14 \ > 1, which 
proves that the quantum correlation thus obtained is non- 
bilocal. 

To study its resistance to noise, one can consider, as be- 
fore, the case where the source sends noisy states of the form 
Pi(vi) — Vi\ipi) + fi)l/4. Because of the non-random 
marginals, the noisy correlation does no longer have the sim- 
ple form of ( 1331 , However, the values of / and J, which only 
involve tripartite correlation terms, are still simply multiplied 
by the global visibility V = v\v 2 . Inequality d23l is thus vio- 
lated for all V > ■ a ■ 

1-fsm Q\ sin 62 

For all values of 61,62, v\ and V2, such that V1V2 < 
1 1 . } — ^-3-, that we tested, we could find numerically an 

l+sin tf\ sin &2 J 

explicit bilocal decomposition for the correlation P^ 4 ^- We 
therefore believe that , , . i — ^-3- is precisely its bilocal vis- 

lH-sm 9i sin 62 ^ J 

ibility threshold Vbu oc , and that once again our inequality d23l 
detects optimally the non-bilocality of f^i . 

Note that as expected, for 6\ = 62 = \ , we recover the 
case of maximally entangled states studied in section El Al 



C. Classical simulation of (noisy) entanglement swapping 



We have shown in Section UlI Al that the correlations ob- 
tained in an entanglement swapping experiment (with a com- 
plete Bell state measurement) can be non-bilocal for visibili- 
ties down to V = V1V2 > 50%. A natural question is whether 
this visibility threshold Vuioc = 50% can be lowered, possi- 
bly by using more measurement settings on Alice and Char- 
lie's sides. 

Studying more complex scenarios, with more settings, 
rapidly becomes very difficult, because of the nonlinearity and 
nonconvexity of the bilocality assumption. We could not find 
so far any scenario where Vuioc could be lowered. However, 
by trying to simulate the noisy entanglement swapping exper- 
iment with an explicit bilocal model, one can obtain a lower 
bound on Vbiioc for all possible scenarios in which Alice and 
Charlie perform projective Von Neumann measurements on 
their qubits, and Bob performs a Bell state measurement. 



1. A fully bilocal model for a visibility V = 25% 



We present here a model that reproduces with a visibility 
V = 25% all correlations obtained in a standard entangle- 
ment swapping experiment, where Alice and Charlie perform 
Von Neumann measurements on their qubits. The model is 
inspired by Werner's model l38ll which reproduces the noisy 
singlet state (so-called Werner state) correlations for visibil- 
ities v = 50%: intuitively, one can simulate the two singlet 
states with visibilities v\ = V2 — 50%, to obtain an overall 
visibility V = V1V2 = 25%; the only nontrivial question is 
how to simulate Bob's Bell state measurement. 

In Werner's model, Alice and Bob share a random vector 
A, uniformly distributed on the Bloch sphere S 2 . After re- 
ception of a measurement setting a, E S 2 , Alice outputs A = 
— sign(a-A); after reception of a measurement setting b £ S 2 , 

Bob outputs B = ±1 with probability p(B\b, A) = 1+B 2 b x ; 
note that this corresponds precisely to the quantum prediction 
for the measurement along the direction b of a qubit in the 
pure state |A) £ S 2 . 

In a similar spirit, we consider the following bilocal model: 
Alice and Bob share a random vector Ai, Bob and Charlie 
share a random vector A2, both uniformly distributed on the 
Bloch sphere S 2 . For measurement settings a, c £ S 2 , Alice 
and Charlie output A = sign(a ■ Ai), C = sign(c ■ A2). As 
for Bob, in order to simulate his measurement, he outputs the 
result B = B°B 1 of a Bell state measurement on a 2-qubit 
pure product state | Ai) | A2), with the probabilities predicted 
by Quantum Mechanics. We show in Appendix ID II that this 
model indeed reproduces the entanglement swapping correla- 
tions, with a visibility V = 25%; of course, one can then also 
simulate smaller visibilities by introducing some additional 
noise. 

We note that there is a significant gap between the upper 
bound V = 50% and the lower bound V = 25% on the bilocal 
visibility threshold, for any choice of measurement settings. 
But the situation is quite similar to the case of locality, where 
Werner's model reproduces the singlet state correlations for 
a visibility v = 50%, while the CHSH inequality Q al- 
lows one to demonstrate nonlocality only for v > 1/ It is 
known however that there exists a local model that reproduces 
the Werner state correlations for a visibility v ~ 65.95% l29ll ; 
it would be interesting to see if that model could be adapted 
to a bilocal model for entanglement swapping correlations. 
The difficulty is to define an adequate simulation of Bob's 
Bell state measurement, with local variables that do not have 
(unlike in Werner's model) a straightforward interpretation as 
quantum states. On the other hand, some inequalities have 
been found that demonstrate the nonlocality of Werner states 
for visibilities v slightly smaller than 1/V2 l39tl . However, 
these inequalities involve a very large number of measurement 
settings. It will certainly be very hard to find better inequali- 
ties to decrease the visibility threshold in the case of bilocality. 
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2. Simulation of entanglement swapping with communication ? 

Toner and Bacon j40ll . followed by Degorre et al. llidll . 
have shown that the use of one single bit of communication 
is enough to classically simulate the quantum correlations ob- 
tained from Von Neumann measurements on a singlet state. 
It is quite natural to wonder whether such a result holds in 
our case, i.e. whether adding some (limited) classical com- 
munication can help to simulate the entanglement swapping 
experiment in a bilocal manner. 

We present in Appendix ID 21 a protocol directly inspired 
from the communication protocol of l4ltl . that uses 2 bits 
of communication. We find that it allows one to simulate 
the entanglement swapping correlations with a visibility V = 
4/9 ~ 44.4%; this is indeed better than the visibility of 25% 
obtained with the previous bilocal model, but this is still a 
pretty low visibility, which does not even reach the threshold 
of 50%. 

It might however be possible to increase this visibility — and 
even obtain a perfect simulation with V = 1 — with a different 
protocol, possibly using more (but still finite) communication, 
in addition to bilocal shared randomness. Recent results in 
this direction indeed look quite promising l42ll . 

V. ON THE ASSUMPTION OF INDEPENDENT SOURCES 
IN STANDARD BELL EXPERIMENTS 

We finally come back to the justification of our bilocality 
assumption. As already emphasized in the introduction, a very 
similar assumption is actually needed in standard Bell tests, 
namely that the sources of randomness used to choose the 
measurement settings are independent from the source emit- 
ting the states that are being measured (the "free choice" or 
"measurement independence" assumption j43l - l45ll ). 

To make this connection more precise, consider a standard 
bipartite Bell experiment, in which a correlation P(a, b\x, y) 
is observed. Bell's local causality assumption writes 

P(a, b\x, y)=J dA p{X) P(a\x, X)P(b\y, A). (48) 

Assume that the random choice of Alice's setting depends on 
the hidden state Ai of her random number generator, so that x 
is chosen with probability P(a;|Ai), and Ai follows the distri- 
bution pi(Ai); similarly, assume that Bob's setting y is cho- 
sen with probability P(2/|A 2 ), where the hidden state A2 of 
his random number generator follows the distribution p 2 (X 2 ). 
The assumption that the settings can be freely chosen implies 
that Ai, A2 and A must be independent. Together with the lo- 
cal causality assumption, the overall probability distribution 
P{a, b, x, y) then writes 

P(a, b, x, y) = P(a, b\x, y)P(x, y) 

= JJJdXiAXdX 2 p 1 (A 1 )p(A)p 2 (A 2 ) 

P(x\Xi)P(a\x, X)P(b\y, X)P(y\X 2 ). (49) 

Let us now compare this situation with the four-partite ex- 
periment depicted on Figure [7] in which Xavier and Alice re- 



ceive particles from a source Si, Alice and Bob receive parti- 
cles from a source S, while Bob and Yolanda receive particles 
from a source S 2 . The four parties perform some fixed (possi- 
bly joint) measurements on their particles, and obtain outputs 
x, a, b and y respectively. In the spirit of our bilocality as- 
sumption, we call trilocal correlations that can be written in 
the form 

P(x, a, b,y)=JJJ dAi dAdA 2 pi(Ai) p(A) p 2 (X 2 ) 

P(x\Xi)P(a\Xi, X)P(b\X, X 2 )P{y\X 2 ). (50) 

Note the similarities with eq. j49l . 

Assume now that for all x, y, the marginal probabilities 
P{x, y) are P(x, y) = P(x)P(y) ^ 0. Under that condition, 
P(x, a, b, y) is trilocal if and only if the conditional probabil- 
ity distribution P(a, b\x, y), where x, y are interpreted as the 
inputs of a bipartite scenario, is local. 

Proof. Suppose that P(x,a,b,y) is trilocal, i.e. that it 
can be decomposed as in 1501 . and that P(x, y) = 
P{x)P{y) ^ for all x, y. From Bayes' rule, we can 
write pi(Ai)P(ar|Ai) = pi(Xi\x)P{x) and p 2 (X 2 )P{y\X 2 ) = 
P2(X 2 \y)P(y). Dividing eq. d50l by P(x)P(y), we find 
that P(a,b\x,y) = P(x, a,b,y)/[P(x)P(y)] is of the 
form d4lfl). with P(a\x, X) = JdAipi(Ai|a;)P(a|Ai, A) and 
P(%,A) = fdX 2 p 2 (X 2 \y)P(b\X,X 2 ) (which constitute 
properly normalised probability distributions). This shows 
that P(a, b\x, y) is local. 

Conversely, suppose that P(a, b\x,y) is local, with a 
decomposition of the form d48b . and that P(x,y) = 
P{x)P(y). Let then Ai and A 2 be copies of the vari- 
ables x and y respectively, so that P{x\X\) = S x \ ± and 
P(y|A 2 ) - 5 y . M . By writing P(x) = £ Ai P(A 1 )P(x|A 1 ) 
and P(y) = p (^2)P(y\X 2 ), and using the fact that 

P(a;|Ai) = S x .\ 1 to replace P(a\x, A) by P(a|Ai,A) (and 
similarly for P(b\y, A)), one gets, from the local distribu- 
tion of P(a, b\x, y), an expression of the form d50l for 
P(x,a,b,y) = P(a,b\x,y)P(x)P(y): hence, P(x,a,b,y) 
is trilocal. □ 

Hence, when Xavier and Yolanda's measurement boxes 
in Figure [7] are interpreted as random number generators, 
which determine Alice and Bob's inputs x and y, the four- 
partite scenario is tantamount to a standard Bell test; the 
trilocality assumption — which naturally extends our bilocal- 
ity assumption — is formally equivalent to the assumption that 
the measurement settings are chosen independently from the 
source S. It is worth stressing that the trilocality assumption 
can thus be tested without any choice of inputs: there is no 
need for any additional "free choice" or "measurement inde- 
pendence" assumption, as it is already explicitly taken into 
account in d50t . 

To finish off, let us illustrate our claims with an explicit ex- 
ample. Consider a four-partite scenario as in Figure [7] where 
the sources Si and S 2 send the (separable) states qi = g 2 = 
i|00}<00] + ||11>(11| to Xavier-Alice and to Bob- Yolanda 
respectively, and where the source S sends singlet states to 
Alice-Bob. 
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FIG. 7. A trilocal scenario, with four parties returning outputs x,a,b,y (here they do not receive any input). Under the assumption that 
P(x, y) = P(x)P(y) > 0, the 4-partite correlation P(x, a, b, y) is trilocal if and only if the corresponding bipartite correlation P(a, b\x, y) 
is local. 



Assume that Xavier and Yolanda both measure X = Y = 
er z , that Alice measures A = |0)(0| <8> o z + |1)(1| ® Ox while 
Bob measures B = g, |0> <0| + ® |l)(l|. The 

correlation obtained by the four parties is 

P(x,a,b,y) = ^[1- , (51) 

leading to 

P( a> b\x,y) = \[l-^=(-ir +b+ ^} (52) 

which is precisely the correlation one would get in a Bell test, 
where the source sends singlet states, where Alice measures 
either a z or <r x , and where Bob measures either CTz +5 x or 

v2 

Jz 7f x ■ This can be understood as follows: the measurement 

of X and Y reveal the inputs x and y of the Bell test; for Alice, 
measuring A = |0)(0| ® <r z + |1)(1| <8> <7x precisely amounts 
to measuring either ct z if a; = 0, or <r x if £ = 1; similarly for 
Bob, measuring B = ^^(g)|0)(0| + ^ z ^ amounts 

to measuring either ° z ^ w if y = 0, or <Tz ^ Tx if y = 1. 

VI. CONCLUSION 

We have developped in this article the study of the bilocality 
assumption in the context of entanglement swapping experi- 
ments. We derived new constraints on bilocal correlations, in 
the form of nonlinear Bell-type inequalities, for different sce- 
narios; in particular, for the experimentally relevant scenario 
where only a partial Bell state measurement can be performed. 
We found in all cases an advantage of the bilocality assump- 
tion compared to Bell's standard local causality assumption, 
as the former lowers the requirements for demonstrating quan- 
tumness in entanglement swapping experiments. 

A lot of questions are left open. For instance, one could 
study scenarios with more inputs and outputs, quantum states 
of higher dimensions, different kinds of measurements per- 
formed by the three parties, etc. It is not easy to develop an 
intuition about which results are to be expected; these ques- 
tions should motivate further work on the study of bilocal cor- 
relations. 

Another natural and very interesting problem is to ex- 
tend our bilocality assumption to more complex topologies 
of quantum networks with independent sources. While it is 



straightforward to formulate the assumption of independent 
sources for the hidden states A (the N -locality assumption) in 
a similar form as eq. (f2|i for instance, translating it into more 
explicit constraints that can be tested numerically, as we did 
in subsection III CI is not trivial; let alone deriving Bell-type 
inequalities for iV-local correlations. For instance, a natu- 
ral extension to our study would be to consider the case of 
an (7Y+1) -partite linear network with N independent sources 
generating singlet states, where all N— 1 parties inside the 
chain would perform some Bell state measurements, and the 
two parties at both ends of the chain can choose among two 
possible projective measurements; based on our numerical 
findings for N = 3, we conjecture that the correlations thus 
obtained can be non- iV-local for overall visibilities larger than 
Vn-Ioc = (1/V2) N (and that this would be obtained by alter- 
nating the horizontal/vertical and the diagonal bases for the 
measurements of each party), although proving it remains an 
open problem. 

Among the various other network topologies worth investi- 
gating in the context of iV-locality, one is especially intrigu- 
ing. Consider a simple triangle with Alice, Bob and Charlie at 
the three vertices (see Figure[8]l. In the quantum scenario each 
of the three edges holds an entangled qubit pair source; the 
three sources are assumed to be independent, hence the global 
quantum state is a product of the form pab ® Pbc ® Pac- In 
the corresponding 3-locality scenario, Alice's output depends 
on the hidden states (Ai, A3), Bob's on (Ai, A2), Charlie's on 
(A2, A3), and the distribution of the three independent states A, 
factorizes. One can thus again easily formulate an adequate 3- 



x 




V 



FIG. 8. A 3-locality scenario, where three parties receive states from 
3 sources, forming a closed loop. Characterizing the (non)-3-locality 
of such a scenario remains an open problem. 
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locality assumption for this scenario in terms of general states 
A,-; however, it is unclear if it is possible to discretize these as 
we did in subsection III CI in order to derive more convenient 
3-locality constraints. Even in the case without any input, 
characterizing the (non)-3-locality in this triangle configura- 
tion seems challenging. 

Studying more deeply the implications of the independent 
sources assumption will lead to a better understanding of the 
nonlocality that quantum networks can exhibit, and of how 
powerful they can be — compared to classical ressources — 
to perform information processing tasks. We also expect 
such studies to lead to new applications, fully exploiting the 
non- TV-locality of Quantum Mechanics. 
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Appendix A: Topology of the bilocal set 

We prove here that although the bilocal set is non-convex, it 
keeps some weaker properties of the local set: it is connected, 
and its restriction to subspaces where the marginal of Alice (or 
Charlie) is fixed is star-convex (which is not the case however 
for the whole bilocal set). These properties can for instance 
be observed on Figures [3] and [4] 

1. Connectedness 

Consider a bilocal correlation P(a,b,c\x,y,z), with a bilo- 
cal decomposition in terms of local variables Ai, A2, and de- 
fine the correlation P^ as follows: each party outputs a result 
according to the probability distribution P (i.e., according to 
P(a\x,Xi), P(b\y,Ai,\2) and P(c\z, A2) resp.) with probabil- 
ity £, and produces a random output with probability 1 — £. 
The transformation from P to P^ is made locally (and inde- 
pendently) by each party; therefore, P^ remains bilocal. For 
£ S [0, 1], P(: follows a continuous path in the bilocal set, from 
P = P$=i to Pq = P^=o, where Pq is the fully random prob- 
ability distribution. All bilocal correlations are thus connected 
to Pq; it follows that the bilocal set is connected. 

2. Weak star-convexity (in certain subspaces) 

Let us consider a subspace of Bp^ a \ x ) C B where Al- 
ice's marginal probability distribution P(a\x) is fixed. Then 
Bp( a \x) is star-convex, meaning there exists a point P G 



Bpi a \x) such that the whole line segment between any P e 

®p{a\x) and P is in B P(a \ x) . 

Proof. Let ii be a product correlation of the form 

P,(a,b,c\x,y,z) = P(a|x)P*(%)P*(c|.z), with P,(a\x) = 
P(a\x), while P+(b\y) and P(c|z) are arbitrary. Clearly, 

P € Bp( a \ x y 

Consider another correlation P 6 Bpt a \ x )> with a bilocal 
decomposition in terms of local variables Ai, A2, and let us 
provide Bob and Charlie with an additional random bit £, such 
that p(£ = 1) = 1 - p(£ = 0) = £ € [0, 1]. Define now the 
correlation Pj as follows: Alice always outputs a result ac- 
cording to P(a\x, Ai); if £ = 1, Bob and Charlie output a 
result according to the probability distributions P(b\y, Ai, A2) 
and P(c\z, A2), resp., while if i = they output a result ac- 
cording to the probability distributions P(fc|i/) and P(c|2i). 
One can easily check that P^ = £,P+ (1 — Clearly, is 
also bilocal, and Alice's marginal satisfies P^(a\x) = P(a\x): 
hence, Pj <G Bp( a \ x ). For £ <G [0, 1], the full line segment be- 
tween P and P in thus in i3p( a | x ), which shows that Bp( a \ x ) 
is star-convex, for the vantage point P. □ 

By symmetry, any restriction $p( c | 2 ) of the bilocal set to a 
subspace where Charlie's marginal probability distribution is 
fixed, is also star-convex. 

The star-convexity property does not however extend to the 
full bilocal set. Consider indeed two bilocal correlations P 
and P, such that both P(a\x) ^ P'(a\x) and P(c\z) ^ 
P'(c\z), and a mixture P ? = £P + (1 - £)P (with again 
£ € [0, 1]). Noting that all bilocal correlations necessarily 
satisfy P(a,c\x,z) = P(a\x)P{c\z), one gets 

P^(a,c\x,z) = £,P(a,c\x,z) + (1 — £)P'(a,c|x,z) 

= £P(a\x)P(c\z) + (1 - §P' {a\x)P' (c\z) 
Pt(a\x)=ZP(a\x) + (l-OP'(a\x) 
P i (c\z)=ZP(c\z) + (l-Z)P'(c\z). 

and therefore 

P^{a,c\x,z) — P^(a\x)P^(c\z) 
= £(1 - Z)[P{a\x) P'(a\x)][P(c\z) - P'(c\z)} . 

For < £ < 1, one thus has P^(a,c\x,z) ^ P^(a\x)P^(c\z), 
and therefore P^ is not bilocal. For all bilocal correlation P, 
there thus exists a bilocal correlation to which it is not con- 
nected by a line segment of bilocal correlations: this proves 
that the whole bilocal set is not star-convex (and therefore, in 
particular, that it is non-convex). 



Appendix B: Explicit local and bilocal decompositions 

We give in Tables Q] to [V] explicit (bi-)local decomposition 
for various correlations analyzed in the main text. Exhibiting 
these decompositions allows us to prove, precisely, that the 
correlation they define is (bi-)local. These explicit decompo- 
sitions are given here in the correlators representation, intro- 
duced in subsection lll C 21 
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TABLE I. Any values of 7 22 and J 22 such that ^T 22 ! + s/W 2 ] < 1 can be obtained by a bilocal correlation P 22 , for instance by the one 
defined by the above explicit decomposition, with K 22 = \/\I 22 — i/j J 22 1. For this decomposition, the constraint \J\l 22 \ + \f\J 22 | < 1 
comes from the non-negativity condition <BU . 

For J 22 = J 22 = \V (and hence K 22 = 0), the table gives a bilocal decomposition for Pq(V) (see section [TTJBTb . valid for V £ [0, §]. 
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TABLE II. Similarly as in the 22-case of Table U any values of I 14 and J 14 such that ^/|J 14 | + ^\J 14 \ < 1 can be obtained by a bilocal 
correlation P 14 , for instance by the one defined by the above explicit decomposition, with K 14 — \/\I 14 \ — \J\J 14 \. Note the strong 
similarities of this decomposition with that of Table [TJ the main difference being that the correlators ej n g with i ^ 11 and k ^ 11 (in the 
fourth sub-table) are no longer internal degrees of freedom of the decomposition; we now display them in non-shaded cells. 
For I 14 = J 14 = \V (and K 14 = 0), the table gives a bilocal decomposition for Pq(V) (section lmATI valid for V € [0, §]. 
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TABLE m. Any values of J 13 and J 13 such that ^|7 13 | + J\J 13 \ < 1 can be obtained by a bilocal correlation P 13 , for instance by the 
one defined by the above explicit decomposition, with K 13 , L and M 13 such that |L 13 +Af 13 | < i(l-i«: 13 ) 2 , |L 13 -M 13 | < 1-(A" 13 ) 2 , 
4|/ 13 | < (l+A' 13 ) 2 and4|J 13 | < i(l-A' 13 ) 2 + L 13 +M 13 < (1-A' 13 ) 2 . 

For J 13 = $V, J 13 = \V, K 13 = ^/ll 13 ] - yiP 3 ! = J% , L 13 = OandM 13 = \V, the table gives a bilocal decomposition for Pq 3 (V) 
(section lllTB~2l >. valid for V G [0, §]. 



For the case where all parties have binary inputs and outputs 
(the 22-case), we use the definition (fl~2l for the correlators. A 
given explicit decomposition will be displayed in the form of 
four sub-tables (see Table HJ, each of them containing the val- 
ues of the correlators e^g, for j = 00, 10, 01 and 1 1; we refer 
to eq. (fl~4i i to clarify the notations on the left column and top 
row of each sub-table. Let us recall that the correlators 
such that i ^ 11, j ^ 11 and k ^ 11 are fixed by the cor- 
relation P 22 we want to reproduce; we display these in white 
cells (the number of which is the dimension of the correla- 
tion space), except for the constant normalization coefficient 
6qoo = 1, shown in a black cell. On the other hand, the corre- 
lators such that i = 11, j = 11 or k = 11 are internal degrees 
of freedom of the different possible decompositions of P 22 ; 
these are displayed in shaded cells. 

Recall also that a valid decomposition must satisfy the non- 



negativity constraint ( fl~6l ): 

for all a, A 7, g(-l) s - ?+ ^- S esjs>0. (Bl) 

This constraint will delimit the domain of validity of our ex- 
plicit decompositions. Finally, the bilocality constraint ([l"7i i 
(i.e. ejQk = e^oo eggs) can easily be checked in each case on 
the first sub-table: the 3x3 bottom-right sub-table (separated 
by dashed lines) must be the product of the column on its left 
with the row above. 

As already observed in subsection III D 31 the 14-case, 
where Bob has one input and four possible outputs, is quite 
similar to the 22-case. The definitions of the correlators, and 
the non-negativity and bilocality constraints on these, are for- 
mally the same. The only difference is that the correlators 
ej u j. with i ^ 11 and k ^ 11 are now accessible experi- 
mentally, and therefore fixed by the correlation P 14 we want 
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TABLE IV. Explicit bilocal decomposition for the correlations Pq 4 (V, rf) of section ITlI C 21 (shown on Figure [5j, for imperfect detection 
efficiencies r/ for Alice and Charlie. We use the notation fj = 1 — r). Three regimes for r\ must be distinguished. For r\ > |, we 
define e v = 0, f v = 1— 77, and the decomposition is valid for all V < V^ iloc = — t-. For | < 77 < |, we define e v — 3~4r], 

/„ = min[4(l-7;) 2 ,4(l-77)(37/-2) + Vr]% and the decomposition is valid for all V < V£ loc = = 4(1-^(2^-1) ^ Finally; for 

?7 < |, we define e v — 2r\~\, f v — Vr) 2 , and the decomposition is valid for all V < V^ iloc = 1. 
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TABLE V. Explicit local (top sub-tables) and bilocal (bottom sub-tables) decompositions for the correlation -Pq 4 (£) introduced in subsec- 
tion IIV AI to study the trade-off between resistance to noise of nonlocality and resistance to noise of non-bilocality. We use the notations 
C{ = cos^^), = sin(£-|), and define = c ^ ^ — - and h% = C£ — S{. 

The local decomposition (which is clearly not bilocal: see its first sub-table) is valid, i.e. dBlt is satisfied, for visibilities 
V < Vioc = l/[c£ + Sj]. The bilocal decomposition is valid for visibilities V < Vuioc = 1/[1 + c^]. 



to reproduce: the explicit decompositions have fewer internal 
degrees of freedom (see TableslHlandllVI-IVli. 

The 13-case, where Bob now has one input and three possi- 
ble outputs, requires slightly different definitions for the cor- 
relators. From a (bi-)local decomposition of a correlation P 13 
in terms of weights q a ^ (with a = aoai, 7 = 7071 = 
00,01, 10 or 11, and /3 = 00,01 or {10 or 11}), we define 
in the 13-case 

e iOk = £a 7 ( — 1)" fc [<7a,00,7 + 9a, 01, 7 + 9a,{10orll},7] ! 
' e ilk = J2a*/(~ I)" fe [<7fi,00,7 + 9a,01,7 ~ 9a, {10 or U} ,7] 1 

, e i2k = Z)a 7 (" 1 ) Q I+, ' fe [9a,oo,7 - 9a,0l, 7 ] , 
which can be inverted into 

fc,oo, 7 = 2- 6 Eis(-l)"' J+ ^ S [e|oS + eilk + ^ 2h } , 

te.oi.7 = 2" 6 £a(-l) a ' W [e M + ej lS - 2e- l2 - k ] , 

'/,,:: ,::, = 2- 6 Ea(-l) a " + ^"[2e^ fc - 2e 7lS ] . 

With these definitions, one can check that again, all correla- 
tors of the form eq^, with i 7^ 11 and k 7^ 11, are fixed by the 



correlation P 3 one wants to reproduce; in particular, one has 
eoo,o,oo = 1 by normalization, and the correlators e- ilk and 
e l,2k w i tn = 10 or 01 correspond precisely to the tripartite 
correlation terms (A X B°C Z ) pi3 and (A X B 1 C Z ) pis b o_ de- 
fined in subsection lIID4l Note that for the decomposition to 
be valid, the weights must be non-negative, which (from 
the previous equations) imposes some constraints on the cor- 
relators ejjk- Finally, the bilocality assumption writes again, 
in terms of correlators, ej j. = ej Q eg j,. 

In Table [TTTJ we display the correlators in three sub- 
tables, using similar conventions as before. 



Appendix C: Comparison of inequality (23) 
with that previously presented in ref. flal 

Here we show that the bilocal inequality previously derived 
in lfl5ll is actually implied by our inequality j23l , for the case 
where Bob has one input and 4 possible outputs. 
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Defining I ± = 2/ 14 ± 2 J 14 , inequality <|23) implies 

I 2 

1+ < 1 + -f ■ (CI) 
Now, from the definitions (|2T1422] |, we have 

j+ = E E 

I-=Y, E PQP^Ewixz) 

with 

P(6°6 1 )E fe o fel (x^) = (-l) fo "^(-l) a+c P 14 (a,6°6 1 ,c|2;,z). 

a.c 

One can see that, up to a small change of notations (bi O 
&o © bi), 1+ corresponds to /, and that |7_ | < 2E, for / and 
E as defined in [15]. Together with ( IClt , we thus find that 
our inequality (l23T l implies the one (/ < 1 + E 2 ) previously 
presented in IU5I1 . 

It is in fact strictly stronger than the previous inequality: 
there exist nonbilocal correlations (such as, e.g., the local cor- 
relation defined by its weights goo,oo.oo = 9oi,oi,oi = h> with 
all its other weights = 0) that violate i neq uality d23l l. but 
which do not violate the inequality of Ref. IU5I1 . nor any of its 
equivalent versions. 



After averaging over Ai and A2, one can easily check that 
all single- and bi-partite correlators in the first 2 lines of ( 1D1I ) 
vanish, as it is the case for entanglement swapping correla- 
tions. To calculate the tripartite correlators, one can show that 
J s2 dAp(A) sign(a ■ A) (A • u) = \a ■ u (with p(A) = i 
being the uniform distribution of A on the sphere S 2 ), so that 
we get 

(AB°C) = J s2 dAip(Ai) sign(a ■ Ai) \\ 

x f s2 dA2p(A2) sign(c • A2) AS. = \ a z c z , 
{AB X C) = \a x c x , {AB Q B l C) = -±a Y c Y , (D2) 

where a X Y Z and c x,Y ' z are the X, Y, Z components of the mea- 
surement settings a and c. 

Now, for the quantum correlations in an entanglement 
swapping experiment, one has, precisely, 

(AB°C)=a z c z , (AB 1 C)=a x c x , 
{AB Q B X C) = -a Y c Y . 

Hence, our bilocal model reproduces the entanglement swap- 
ping correlations with a visibility V = j = 25%. 



2. Simulation with 2 bits of communication 



Appendix D: Classical simulation 
of noisy entanglement swapping 

We give in this appendix the details of the simulation of 
noisy entanglement swapping correlations by bilocal models, 
as presented in section llV CI 

1. Simulation without communication 

Consider the bilocal simulation protocol presented in sec- 
tion llVCll To study the correlation P(A, B°B 1 ,C\a, c) it 
defines, it is sufficient to calculate the average values 

(A),(C) > (B°) > (B 1 ),(B°B'), 
{AC), (AB°), (AB 1 ), (AB^ 1 ), (B°C), (^C), {B a B l C), 
(AB Q C) 7 (AB 1 C) 7 (AB°B 1 C) . (Dl) 

For B = B°B 1 e {++, H — , — h, } being the result 

of a Bell state measurement, corresponding to the outcomes 
|$ + ), \ ^ + ) and \^f~) respectively, it is convenient to 
note that the outcome B° actually corresponds to the mea- 
surement of <7 Z £3> <7 Z , the outcome B 1 corresponds to the mea- 
surement of a x <S> <7x, and the product B°B 1 corresponds to 
the measurement of — <j y <8> <r Y ; therefore, the expectation val- 
ues for a given quantum state | Ai) | A2) are (B )^ r = \\ A 2 ,, 
(B 1 )^ = A X A X and (B^)^ = -A Y A Y (A*' Y ' Z being 
the X, Y, Z components of Ai). 



We have mentioned in Section IIVC2I that one could in- 
crease the visibility of the simulation with the help of commu- 
nication. Inspired from the communication protocol presented 
in 14 111 (Theorem 10) that simulates the singlet state correla- 
tions, we slightly modify the previous bilocal model in the 
following way: instead of starting with Ai and A2 uniformly 
distributed on the Bloch sphere, we use 1 bit from Alice to 
Bob, and 1 bit from Charlie to Bob to bias the distributions of 
Ai and A2 according to (see Theorem 6 of l4lll ) 

Pa(Ai) = -^-, P&*) = -^ r . 

One can still easily check that the single- and bi-partite cor- 
relators in dDU vanish, as for the entanglement swapping cor- 
relations. The tripartite correlators can be calculated in a sim- 
ilar way as in (ID2t . with the modified distribution functions 
Pa(Ai) and ps(^2) above. An easy calculation shows that 
X52 dAp a (A)sigii(a ■ A) (A • u) = ^ f £2 dA(a ■ A) (A • u) = 
I a • u, from which we now get 

(AB O) = la z c z , {AB l C) = § a x c x , 
(AB°B 1 C) = -|a Y c Y . 

This shows that our bilocal model, augmented by two bits of 
communication, reproduces the entanglement swapping cor- 
relations with a visibility V = § ^ 44.4%. 
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